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THE gq? EQUATIONS WHOSE SOLUTIONS ARE 
PRODUCTS OF SOLUTIONS OF q¢’ EQUATIONS 
OF LOWER ORDER 


By F. H. JACKSON (Eastbourne) 


[Received 4 July 1939] 
(1) 1. Introduction 


Ir is well known that, if we have two homogeneous linear differential 
equations of the second order 


y" +filx)y’ +fo(x)y = 9, 
n” +(x )n’ +¢44(x )n = 0, 


then a linear differential equation of the third order can be found 
which will be satisfied by the product yx 7, conditional upon an 
equality between the invariants J, formed from the coefficients f, ¢, 
of the two equations of the second order. 

In this paper certain general theorems for q-product difference- 
equations will be found. These theorems are the g-function analogues 
of the above theorem, they are obtained by very simple analysis, 
and contain the theorem referred to above as a particular case. They 
may be worth placing on record. Examples will be given dealing 
with the products of g-hypergeometric, and other functions of various 
orders. In the case of the g-hypergeometric equation of the second 
order twenty-four integrals will be obtained corresponding to the 
twenty-four Kummer integrals of the ordinary hypergeometric dif- 
ferential equation. The paper is divided into two parts. Part I is 
devoted to analysis of four general forms of g-difference equations. 
Part II contains special examples. 


2. Notation 
[ use a, to denote a functional coefficient a(x), and, in general, 
a,, to denote a(qg"2). Then 


{q¥+-a)q°+b jy = 0 


is the general g? equation of the second order, analogous to the 
homogeneous linear ordinary differential equation of the second 


order. Here q® is an operator exp(log a2 7) which, when operating 
x 


3695.11 








2 F. H. JACKSON 
on any function ¢(x) (possessing derivatives) converts ¢(x) into 4(qz), 
as is easily seen if g? be expanded in an exponential series and the 
infinite sequence of operations denoted by 0, 6?, 6°,... be carried out. 
Throughout the paper factors such as g°+*—1, q?8+*—1 will be 
denoted by [@+<], [(20-++a)]. Special g? equations such as 
{ 26 | ax+b a 


lq- C = 
\? cx+d- ca+ 


can be transformed into 
Aa[O+o][O+B]y—[]O+y—1ly = 0, 

very conveniently for analysis. 

The equation {g79 1 ag? +-b+cxy = 0 

jue qq rag wsy ’ 
in which a, b, c¢ are constants, contains implicitly Bessel’s equation 
and others allied to Bessel’s equation, and also others which might 
seem to have no connexion with Bessel’s functions; for example, 
2 O 1 2 : 
(q??—q?+2)y = 0 

has solutions in series first discussed by Rogers and later by Rama- 
nujan. This equation is satisfied by 


2 4 
_ a ia. a a 


gig—1)  g(¢?—1)(q*—1) 


which by inversion of g becomes 
2g VA 

a vG 

ce at a 

1—q* = (1—q*)(1—q") 

3. THeorEM A. A q equation of the fourth order can be found 
which is satisfied by the product of the respective solutions of two equa- 
tions of the second order, conditional on an equality between certain 


12 


ae 


invariants of the two equations of the second order. 


Let 


{q*?+-a, q' +bo}y = 0 | 
{920 eS (x) 
{q"’+c99°+d}n = J 


be equations satisfied by y = ¢(x) and » = (x) respectively. For 
the sake of brevity we shall write 


o(q'x) = ¢(r) (r = 0,1, 2, 3). 
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We obtain from (a) the following six ge 

$(2) +49 $(1)+-by (0) = (1) 
$(3)+-4, $(2)+6, (1) = 0, (2) 
$(4)+4, $(3)+6, ¢(2) = 9, (3) 
(2) +-Cy (1) +d H(0) = 0, (4) 
(3) +c, p(2)+-d, (1) = 0, (5) 
(4) +, (3)+-d, (2) = 0. (6) 


On combining equations (1) and (4), and placing all oddly-paired 
terms on the right-hand side of the resulting equation, we obtain 


$(2)yb(2) +g Cy H(1)yh(1)—by do (O)yb(0) = —ey (2)9h(1) a9 (1)yh(2). 
(7) 
Similarly, from (3) and (6) we form 
(4 )yb(4) —ay Cy 6(3)xb(3) —b, dy J(2)xh(2) = by cy 6(2)h(3) +a, dy $(3)xh(2). 
(8) 
Pci (7) by an arbitrary coefficient A, and then subtracting 
7) from (8), we obtain an equation in which all the terms on the 


oft hand side are evenly paired, but on the right the expression is 


’ wis Ay 1)\ ovbareys 40 5 
baeabl2)H8) + 52 HUD] anda (2), 6(3)+ 2°9- #0) 


This expression, by suitable choice of A, on comparison with equa- 
tions (2) and (5), is convertible into 
— {be Cy Coty Ay dy}h(2)xp(2) 
provided that A = byc,d,/cg = da, b/d, 
so that we obtain finally 
(4 )yb(4) —Ay Cy 6(3)eh(3) —b, dy f(2)xh(2 Sil )eb(2) —Ady Cy H(1)y(1) + 
+Abo do $(0 on = —{¢,C,b. +a, 4, d,}h(2)yp(2) 
This, when expressed in terms of g°, fa is 
(q29 —ayc, 29+ (Cy Cy bg +A, Ay dp—by dy—A)q*? —Nay Coq? +-Abyd JY = 0 
(9) 
a linear equation of the fourth order satisfied by 
Y = $(x) x p(x) 
a product of solutions of the two equations of the second order. The 
condition stated above may be written 
I = ayb,/ayb, = c,d,/¢g dy. (10) 
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These functions of the coefficients are invariant for the following 
change of the dependent variable in equations of the second order. 
On replacing y by 7 Xf in the equation 

{q9+-a,q°+boly = 0, 
we obtain Nefot fi +59 = 0, 


Sie 5," So) 


which is q?+a3q°+ y =: @, 
"fe ho) 


an equation of the form 
204 ta 
{(g+Aog’+ Bo}n = 0, 

in which 

Ay = fifo, A, = %f./fs, A, = a2fs/f,; 
and similarly, 

By = bofolfe, B, = by filfs; By = be folfa; 
from which it follows that we have the invariancy 

A, B,/A, B, = a,6,/a,6, = I. (11) 

As in the theory of ordinary differential equations, this invariant 
property is of great utility, enabling us to transform equations into 
others more amenable to analysis, especially in cases where direct 
solutions of the original equation fail, owing to divergence of series, 
or are impracticahle because of difficult recurrence relations between 
successive coefficients in the series. Examples illustrating both these 
cases are shown in the second part of the paper. 

4. The question arises: ‘May not the product of solutions of the 
two equations of the second order satisfy an equation of the third 
order, simpler than that of the fourth order?’ The answer is in the 
affirmative, as is easily seen by combining 

$(2)-+ay $(1) +b96(0) = 0 
with b(2)+-c¢y (1)+d, (0) = 0, 
and $(3) +a, $(2) +5; (1) 
with (3)+-c, (2)+-d, f(1) 
Then, using an arbitrary coefficient A to eliminate odd pairs, we 
obtain 


{a"?— a (ay a,—b,)q??-+-d, (ay a, —b,)q°—Aby a, Y=0 (12) 
0 


satisfied by Y = d(x) x (2) 
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subject to the condition 


A = a, d,/Cy = by ¢4/ap. 


Inspection of this condition shows us that, if the product of solutions 
of two equations of the second order is to be such as will satisfy an 
equation of the third order, then a closer degree of symmetry must 
exist between the coefficients (@», by), (cy,d,) than is necessary if an 
equation of the fourth order is to be satisfied. For example, in the 
case of the qg-product of two basic Bessel functions J,,(q, x) .J,,(q, x) 
such product satisfies an equation of the fourth order, but, ifm = n, 
then the simpler equation of the third order will be satisfied.* 
We may note that, when a,a, = b,, then cyc, = d,, and in such 
case equation (12) reduces to 
{q—aya,4_,Cy¢,¢_}Y = 0, (13) 
satisfied by the product of solutions of 
{q?+-a9q°+a,a-,}y = 0, 
{q9+-¢99°+¢ye_s}n = 0, 
from which some results of interest follow. I confine this paper to 
the equation of the fourth order as being more general and as con- 
taining implicitly that of the third order. 
5. THEoreM B. If 
{q+ q%+by}y = 0 
be satisfied by y = ¢(2), 
then will {q*9+-a, g°—ayb, q?—by bY = 0 
be satisfied by Y = $(x) x $(q2). 
Using the notation of §2, we have 
$(4)+-ay $(2)+-by 6(0) = 9, 
$(5) +a, $(3) +b, (1) = 0, 
whence 
$(5)$(4) +a, $(4)6(3) = ay byG(2)4(1) +55, 4(1)4(0), 
which is {q*9-+-a, g?9—ayb, g? —by b,}h(qx)o(x) = 0. 
This theorem, obtained by most elementary arithmetic, depending 
on a property peculiar to the first six integers, is remarkable as con- 
taining a great number of interesting theorems regarding products 
of hypergeometric functions. Such products have been discussed by 
* F. H. Jackson, Proc. London Math. Soc. (2), 2 (1904), 207. 
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several writers* in the case of F(a,8;y;x). The g-theorem (B) was 
given by me some years ago. I return to it here as falling into line 
with the rest of the paper, also because of the simpler notation, 
and because of an example much more general than those previously 
discussed. 

6. THroREM C. This theorem, which is concerned with the general 
q° linear equation of order n, is useful for transforming such equa- 
tions into forms capable of solution in hypergeometric series. We 
consider first the simple case of an equation of the second order, 

f 20 6 ; 
{q°+ay9q°+bojy = 0, 
satisfied by 4(x), and we associate with this an equation of the first 
order, 6,6 2 
iT —Cosn = 9, 
supposed to be satisfied by 
ob(2). 
Then, since 
(qx) = cyrb(x), 


p(q'x) = ¢, p(qe), 


b(q"x) ‘ Cy ¥(q" 2), 


we see that {q¥+-ayc, gq? +bycycs}¥ = 0, 
must be satisfied by 
Y = yX7 = $(e) x(a), 
and it is obvious for an equation of the nth order, by similar reason- 
ing, that, if an equation 
{q+ Aggy” + Bog” +... + Koly = 0 
be satisfied by y = $(2), 
then, on associating this with an auxiliary equation 
{q°—co}n = 0, 

we can form 

{q?9 4 Ay Cy _1 "Mt Bg Cyn—a Cn—2 FM ... +09 C}---Cn—1 Koh Y = 0, 
which must be satisfied by 

Y =X». (16) 
* Watson, Quart. J. of Math. (1907); Appell, Comptes Rendus, tom. 91; 


Goursat, Annales Sci. de ’ Ecole Normale, Ser. 2, tom. 12; Jackson, Messenger 
of Math. 474 (1910), 92; Nicholson, Phil. Mag. (1905). 
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3y suitable choice of cy in the auxiliary equation it is often possible 
to reduce equations of type (15) into g-hypergeometric forms. 

(II) 7. This part of the paper is concerned with special examples 
of the theorems obtained in part I. The general equation of the 


fourth order is 
(q#9 —ay Cy G*° + (Cy Co bg, Ay dg—b dg—A)q”?—Ady Cy q?+-Abydy}¥ = 0, 
| (17) 
satisfied by the product of solutions of the respective equations 
(q*? +a q° +boly = 0 
(q°+-coq’+do}n = 0, 
in which dp, bo, Co, dy satisfy the invariant equations 
I = a,b,/ayb, = cyd,/cyd, 
and A = bycyd,/cy = dyayb,/ay. 
If now we choose a, = A, b, = B(1—2?*) and A and B are constants, 
then, in order that the invariant equations may be satisfied, we see 
that we have two choices for Cy, do: either cy = a, dy = by, or 
alternatively 
Cy = C/(1—q?x?), d, = D/(1—q?z?). 
We take the latter case only, so that we have 
b, = B(1—z?), b, = B(1— q*2"), b, = B(l—¢q*z*), 
Cy = C/(1—q?x?), ¢, = C/(1—q'*z?), Cc, = C/(1—q®z*), 
d, = D/(1—q?x?), d, = D/(1—gq'*z?), d, = D/(1—q*x?). 
On substituting these values in the general equation (17) we obtain 
a2{qi9+6___ BD)(q?9+4+4-9?29+2)4 B2D yy — 
—{(q49—2 BDq**+ B2D?)—q*(Aq®— BC)(Cq®—AD)}¥Y = 0, 


which, by a well-known device, can be reduced to q-factorial form, 
on changing the constants A, B, C, D, as follows: 


A=q"+q", B=", 
C= g+q’, D = qe’; 
and at once the second-order equations take the forms 
[6—m][6—n]y—a*y = 0, 
qh +[O—p][O—v]n—2*g*?+*n = 0. 
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The fourth-order equation takes the form 

x*{ (20-+-4—m—n—p—v)|[ (20+ 2—m—p—n—v)|¥ — 

—[06—m—p][0—n—p][O0—m—v][06—n—v]¥Y = 0. (19) 

The doubled brackets are convenient as showing us that in solving 
by power series the factors arising from such terms of the equations 
will advance by q?. The doubled brackets will be used in the series 
wherever advance is by q?. 

As the method of solving these equations is similar to that used 
in solving ordinary equations, it is unnecessary to enter into detail. 
Inspection of the three equations shows us that we shall have indicial 


equations [p—mJp—n] = 0, 


[p—r]lp—v] = 0, 
[p—m—p]=0, [p—n—p] = 0 
[p—m—v] = 0, [p—n—v] = 0, 
if solutions of the form >} A, x"*? are sought, and thence a funda- 
mental solution for the equation of the fourth order will be, with 
four arbitrary constants, 
Y = Ur" +¥O{[ (4+-m+p—n—v)]| (2+-m+p—n—y)]; 
[2+-m—n][1+p—v][2+m—n+p—yv]; 2?}+ 
+Ba"t+O,+Cx™+'O,+Da"+"D,, (20) 
where ®,, ®,, ®, are obtained from ® by interchange of (i) m and n, 
(ii) » and », (iii) m and n, p and v respectively. 
The solutions of the equations of the second order are 


\ 4 pan! a | 
$e} ba" 1+ 05 atop 


— mf ya 
note § (2]2+n— 


Bystacal 


| x*g?te- v ) vf} la ) 
” = ce 1+ +...}-+D2” 
( [2] 2+u—-r] "J "BRP 
The solitary q-factors in the numerators of the series y are, for the 
(r+1)th terms, 


1+ 


ghu—)+er. gre-H+2, 
If, for brevity, we denote the series by their initial letters we can 
form the following identities: 

(W)=(a.c), (B)=(b.c), (©)=(a.d), (D) = (b.d). 
Of course the arbitrary constants must be suitably chosen and will 
usually be unity. 
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A similar set of four identities might have been obtained from 
solutions of the equations in series of descending powers of x. Where 
such solutions fail owing to the presence, after a finite number of 
terms, of zero factors in the denominators of terms, solutions can 
be found by the substitution of series of the form 


> {A n x” +} B,, x” +P log x}, 


[ omit discussion of such solutions. The object of the present paper 
is to introduce the general theorems given above and to show a con- 
venient notation for forming a general theory of q° ‘differential’ 
equations. 

8. An interesting example of the general theorem (B) is afforded 
by certain g-hypergeometric equations of the fourth order. The case 
of the ordinary hypergeometric series has been discussed by several 
writers (loc. cit.). The series then discussed had two independent 
elements only, since in F'(2a, 28; 2y;2) the condition 2y = 2a+28+1 
was imposed. In the following example the series will have four 
elements subject to the condition 


28 = 2a+2B—2y—1. 


The general theorem is that 


{q*°+-a, q*°—ayb, ¢?—by bY = 0 (20) 
will be satisfied by Y = d(gx)d(2), 
where ¢(2) is a solution of 
(qa q%-+bo}y = 0. (21) 


Consider the special equation 
afag*+1+ bg? + Lyy—{aq*®+-cq??+ ly = 0, 

in which a, b, c are independent constants. This equation will have 
solutions of g-hypergeometric form, for on replacing the expressions 
in the brackets by pairs of factors, the hypergeometric form is 
obvious. We see also that on comparison with the general form (21) 

ay = (bx—c)/(aqgr—a), by = (w—1)/(aqx—a), 

a, = (bgx—c)|(ag2x—a), b= (qe—1)/(ag*e—a). 
On substituting these values in 

(q’+a, q*? —ayb, P—by bY =0 

we obtain 


a(ag??+1— 1)(ag*®?+1+-bg?+1)¥ —(ag®—1)(ag?-+-cq?+1)¥ = 0. 
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For the sake of symmetry we change the constants a, b, c as 
follows: 

Replace a by g20+?8-1 — g2v+28, 

b by —q?*—q-F, 
Cc by —q?”7—q-* 
with the condition 26 = 2a+28—2y—1, leaving the following equa- 
tions, with three independent elements: 
al (20-+2«)][ (20+ 28) y—[ (20+: 
gx{ (20+ 2a) ][ (20+ 28) ]}n—[(20-+2 
a| (20+ 20+ 28) |[6+ 2a][ 0+ 28|Y — 
—[(204- 274-28) [0+ 2y][6+28]¥ = 0. 

The fundamental solutions of these equations in g-hypergeometric 
power series can be written down on inspection. For series of the 
form > A, x04 : 
we have indicial equations 

[2p+2y] = 0, [2p-+28] = 0; 
[2p+2y+26]=0, [pt+2y]= 9,  [p+28] = 0, 
the principal roots being 
—Y, —$; —y—6, — 2y, — 26, 
respectively, so that, denoting arbitrary constants by A, B, C,..., we 
write (0 denoting g-hypergeometric series) 
y = Ax-7{[20—2y][28—2y], [28—2y+2];2}+ 
+ Bx*O{[ 2a— 28][ 28 — 28], [2y—28—2]; 2}, 
1 = Cx-¥{[20—2y][28—2y], [23—2y +2]; gzx}-+ 
+ Da~*O{[ 2a— 28 ][ 28— 28], [2y—28—2]; gz}, 
Y = Ex-?O{[ 20+ 2B —2y— 28](2a—y—8)(2B—y—S), 
(y—8+1)(8—y+1);2}+ 
+ Fa-*{[2a-+ 28—4y](2a—2y)(2B—2y), 
[28—2y+2](25—2y+1); x}+ 
+ Ga-®@{[ 2a-+ 2B — 48](2a—28)(2B—28), 
[2y—28+2](2y—28+41); a}. (22) 
The square-bracket factors (in the series for ®) advance by q?; the 
factors in round brackets advance by q. 
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We can equate identities from products of pairs of the series for 
y and 7, with single series taken from Y, namely 


(AxC)=(F), (BxD)=(@), 

(Bx C)+(AxD) = (£), (23) 
provided that the arbitrary constants be suitably chosen (usually 
unity). Three similar theorems for the fundamental solutions in 
descending powers of x can easily be written down. I omit these 
for the sake of brevity. The theorem is more general than that given 
previously by me,* which contained only two disposable elements 
x, 8 as against three in the present example, namely, a, B, y, 6 
subject to 26 = 2a+28—2y—1. 

9. In this section of the paper some examples of Theorem C are 
given, which will illustrate the use of auxiliary equations. 

The following auxiliary equations will be required in subsequent 
work. They are given here as lemmas. 

Lemma I. The equation 

{q—9’x"}n = 0, 
in which p and c are constants, is satisfied by a function 
7 = S(p,c, 2), 


1 ie (1+ q”)? ® 
wh "e S c.Zt) = a . " ° (25) 
ere (p Cc r) (1+ gPa*) | (1+q”a%q?)(1+q”"a—q-?) 


For, on replacing x by gx in this function, the function repeats with 
a multiplier g?a*. In various notations this function is well known. 
Jacobi gave its expansion in a Laurent series. It is closely allied to 
Weierstrass’s sigma function, as is obvious if x be replaced by exp é. 
[ have discussed it previously in connexion with the basic gamma 


and double gamma functions, etc. 
Lemma II. The equation 
{q’—(x—a)}y = 0 
has the solution 
x s 
q Popes - 
7 eT ata) Qa) 


. Jackson, Messenger of Math. (1910), 92. 
. Jackson, Proc. Royal Soc. A, 76 (1904), 127. 


rhe 
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For, if we substitute } A, «+e for y, the recurrence relation is 


found to be Ani(q"tei+a) = A,, 


and the indicial equation is g°+a = 0. 

Similarly, {q?—(x—a)-}}n = 0 
can be solved. These series are well known. 

10. To solve {9794 rx-4q9 + ua }y = 0. (27) 

Direct solution of this equation by substitution of Laurent series 
for y is either impracticable or impossible except in such simple cases 
as a/b = +1, 2, or $. When the ratio of a to b is any number other 
than the above, the recurrence relations are such as to make any 
calculation beyond one or two terms impracticable. When a/b is 
unspecified numerically, it would appear impossible to obtain recur- 
rence relations between successive coefficients. In the very simple 
case A = —p = 1= a = BJ, solutions in the form of Rogers’ series 


. gra” 
[~) ” 
are easily obtained. 
To obtain a general solution of (27) we use an auxiliary equation 
{q?—2"}n = 0. (28) 
On comparing (27) and (28) with the general forms 
(q*+ayq’+bjy=0,  {9’—c}n = 0 


we see that P 
' {q@+ayc q+, Co ey}yn = 0, 


Co = 2, Gi, = Az, 
*£= fx @, = Ag *s*, 
so that {q?9 + Aae—4q9 +e 4 gra%—Yyn = 0. 
Two cases arise, and both can easily be solved in Maclaurin series: 
(i) c = b—a, (ii) c= a. 
Case (i) gives us 
{q?9+ rae—ag9*e+ gear yn = 0. (29) 
Substitute $ A, 2?-"¢-® for yn. This gives an indicial equation 
g? = —p/A, and a recurrence relation 


A s 419 (Get 1) —_ peg?na—c) 4 A, 
from which, on writing 
é for wat-“[Aty?,  Q for g 











ON @ EQUATIONS 
we have a solution 


mses SEG) 


The auxiliary equation 
{q’—2*}y = 0 
by Lemma I is satisfied by 
— S(q, x), 
so that, finally, 
Qr’- =e 
(Qr— 
The series and product are both a if |Q| <1. 
Case (ii) gives us the equation 
{q°?-+-Ag?+#+ pga hyn = 0. (31) 
Solving by substitution of } A, 2"@¢—+ for yn (with the auxiliary 
equation {q’—2z%}n = 0), we obtain 


y= Ag2*{1+ re = , Q-min 2) — S(q, x), (32) 


+S(q, x). 


y=A oa {+ 


in which, as before, 
ge = —)Aq* Q — g-, E = puac?a—b /d2q4, 
The series and product are convergent if |Q| > 1. 
As previously remarked, these solutions reduce to series of the 
form of Rogers’ series, when b = a = A = —p = 1. 
It is, I think, worth noting that in the special case 


y= (33) 
direct solution by substitution of a Laurent series for y gives recur- 
rence relations connecting sets of three coefficients and a solution 
y = {1l+axr+q(a—1)a*}—{q3x3 + q®axt+q(a—1)a5}+-... ad inf., 

—{(a—1)a-!+-agq?x-*+-g*x-4} + {(a—1)q'x-4+- gaa + ga} — 
—...adinf., (34) 
in which all the indices of qg are triangular numbers, and a is the 


ratio A, Ay, of the first two positive coefficients in the Laurent series 


+0 
> ae 


n=— © 
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substituted for y: the general terms are easily written down. An 
investigation of the properties of this function, a quasi-0-function 
might prove of interest. I have not met with any series of this form 
previously. Some interesting cases arise by giving special values to 
the arbitrary constant a. 


11. Equations such as the following are easily solved: 


{q?94 (ax+b)q9+ (ca?+da+e)hy = (35) 

{929+ (ax+b)-199+ (ca?+da+e)y = 0, (36) 

{q*9+-aq?4 + b—a*ly = 0, (37) 

{(Aat+ Bab)qn9+ (Cat-+Dat)g-+....+(Jat+Kak)}y = 0 (38) 

in which the indices (a, ¢, e,..., i), (b, d, f,..., k) are numbers belonging 


to one and the same arithmetical progression. It will suffice to 
indicate by one example how solutions are obtainable through use 
of an auxiliary equation. 

Write equation (35) in the form 

(q?-+ (ax-+b)q9+-c(x—p,)(x—ps)}y = 0, 
and form an auxiliary equation 
{¢—1/(e—p,)}n = 

then from the general equation we obtain 

{qg294 (aax-+-b)/(qx—p,)q?+e(a—p,)(w—ps)/(qu—p,)(w—p,)}Y¥ = 0, 
ax+b ? ¢ the \y 

=~," ° qe —piJ 

This equation (since the coefficients have a common denominator) 
can immediately be solved in series of the g-hypergeometric form. 

The » equation is solved (Lemma II), hence a solution y, = ¥,/n, 


reducing to i+ = 


is easily found. 
By using another auxiliary eae 
{¢’—1/(x—p2)}n = 0 
a second particular solution of ip is obtainable. 

It is unnecessary to write out this solution which merely inter- 
changes p, and p,. To solve equation (36) the auxiliary equations 
will have to be of the type 

{q’—(x—p))}n = 0, — {g’—(a—p)}n = 0. 
Solutions of these follow from Lemma = so eis two solutions of 
(36) in the form of g-hypergeometric series divided by y follow at 
once. 




















ON gq? EQUATIONS 
Equation (37) requires an auxiliary equation 
{q—e/(1—x)}y = 0, 
which will give 
{q*9 + acq?4 + be} ¥ —a{g*?+2#—be}¥ = 0, 
on putting 
a= —q-"—q-*", b = q-m-2n, c=q™. 


Solutions analogous to 


io * { - (2n+-3) 25 ) 
= mpi) ) * 2@n+2) Tf 


are obtainable. Also, by changing q? to g and putting 


exp(ix)J,() 


n—m = 2p-+1, e = CB, 


we get 


In _lL 252 ~ 2,2 
of a = ee i ' —a-at} x ,F{[2p+2], (2p+1); cz}, 
where round brackets advance by q, square brackets by q?. 

Equations of type (38) can be solved by means of an auxiliary 
equation {q?—2-jn =0 (A=c—a = d—), 


giving an equation of the form 


n 


n 
> Y,. gea—2grO > %,. sty) F 
r=0 r=0 


capable of solution in g-hypergeometric series in powers of 24", since 


the equation can be written 


(x 5 a, g? + 3%, q°}Y = 0. 
The particular solutions will be of the type 

nr ,{%4; Casesey Map B,, B, : g-"), 
There will be such particular solutions from the » principal roots 
of the indicial equation = 


> Bp" = 0. 
r=0 


It would be tedious to write out at length such solutions, which 
the reader can obtain for himself on lines parallel to the case of 
solutions of ordinary differential equations of similar type. The 
object of this paper is not to give a complete theory of such equa- 
tions, which would demand a treatise, but rather to introduce a 
notation and method for obtaining solutions of large classes of q? 


equations. 
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12. The twenty-four integrals of the g-hypergeometric 
equation of the second order 





The equation is 
nq{0-+o](0+2ly—[6l0-+y—1y = 0 = y—a—P). 


In the following work I write 





=| Ti (1—2q") 

: ] (1—ag"*)’ 

The integrals are analogous to those given by Kummer for the 
ordinary hypergeometric differential equation. As these integrals 
have not been given previously, except in four cases,* I propose to 
state fourteen of them here and to indicate how the remaining ten 
integrals can be written down: the numbers will follow Kummer’s 








order. fa, B; y3y2/a8)], (I) 
with a transformation of this g-hypergeometric series into 

(1—2), ®[y—a,y—B;y;2]... (Heine), (IIT) 

(1—x)_, O[a, y—B; y; gu(u—q**)] _ (XVII) 


\ 
} 
(1—2x)_,®[B, y—a; y; gu(x—q*1)] | (Jackson). (XVIII) 


Fe ofa, Bia+B—y+1;q(1—2)) (Vv) 
T(y—a)LG(y—B) 
Transformation (V) is only a special case of a general theorem pre- 
viously given by me.+ 

As in the case of the ordinary differential equation, inspection of 
the gq’ equation enables us to write down the following particular 
integrals in positive or negative powers of 2: 


a1-¥O[a—y+1, B—y+1; 2—y; ya/aB], (11) 

2-*D[a, a—y-+1; a—B+1; 9/2], (IX) 

a$@[B, B—y+1;B—a+1;q/z]. (X) 
Transforming (II) by means of (III), we have 


x1-7(1—z), O{1—a, 1—B; 2—y; a}. (IV) 


* Heine, Kugelfunktionen (1878), 99. F. H. Jackson, Proc. London Math. 
Soc. (2), 2 (1904), 209; Messenger of Math. (1908), 62; Quart. J. of Math. 
(1927), 377. 

{+ F. H. Jackson, Messenger of Math. (1908), 62. 
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Transform (II) by (V), then 
xt-(1—a), ®fa—y+1, B—y+1; a+ B—y—1; g(1—aq’)}. (VI) 
Transform (IIT) by (V), then 
(1—2), ®{y—a, y—B; y—a—B + 1; g(1—2q")}. (VII) 
Transform (IV) by (V), then 
ax1-v(1—2), O{1 —a, 1—B; y—a—B +1; (1—aq)}. (VIII) 
Transform (II) by (XVII), then 
x1-¥(1—2),_4-1 P{a—y+1; 1—B, 2—y; qu(w—q>-7**)}. (XTX) 
Interchanging a, 8 in (XIX) we obtain (XX) 
The remaining ten integrals can be obtained by transforming the 
negative power series (IX) and (X), using transformations (III), 
(V), (XVII), (XVIII), and noting that ga-! interchanges with 
q*+8-v+1z-1; also that the factors (1—aq*) advance by q in each term 
of the series, and are not powers of (1—2q*), where A = y—a—fB. 








EXPANSIONS AND INTEGRAL-TRANSFORMS 
FOR PRODUCTS OF LAGUERRE AND 
HERMITE POLYNOMIALS 
By E. FELDHEIM (Budapest) 

[Received 17 July 1939] 


THE present note contains various problems concerning Laguerre 
and Hermite polynomials. I begin with an expansion of a product 
of Laguerre polynomials in a series of Laguerre polynomials and 
establish the expansion of a product of three Hermite polynomials 
in a series of Hermite polynomials. Then I construct the Gauss 
transform of a product of two Hermite polynomials and give some 
formulae in connexion with this transform. 

1. I shall consider in the first part of this note a problem similar 
to that discussed in a previous note:* the expansion of a product 
of two Laguerre’s polynomials in a series of such polynomials. In 
the note mentioned above, I obtained the developmentt 


L(x) L(x) = + a (—)sBOrr LO(r) (m>n), (1.01) 


m 
s=m—-n 


where the generalized Laguerre polynomial is defined by the formulae 


; ; ya Xer d®(e-*y" | x) n Nn + Y (—a)™ 
ealeaie n! da” — 2 n—m) m! ~ ae 
m=0 


The coefficients BY”"™ in (1.01) are calculated in the following 
manner: we expand the product on the left of (1.01) in a series of 
powers of « by means of (1.02) and, replacing the powers of x by 
the corresponding expressions given by 


— s—™m 
m=0 


we obtain finally (1.01). For the explicit expression of BO", see 


my note.* A. Erdélyit also established by the same method a 


* KE. Feldheim, ‘Quelques nouvelles relations pour les polynomes d’Her- 
mite’: J. of London Math. Soc. 13 (1938), 22-9. 

t Cf. also G. N. Watson, ‘A Note on the polynomials of Hermite and 
Laguerre’: J. of London Math. Soc. 13 (1938), 29-32. Watson’s method is 
based on the generating function of Laguerre polynomials. 

{ A. Erdélyi, ‘On some expansions in Laguerre polynomials’: J. of London 
Math. Soc. 13 (1938), 154-6. 
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generalization of our problem, giving for the coefficients of the 
expansions 

Li) (ley 2) LS (Iby 2)... LS%2)(Iey 2) 


M+ M,+ ...+My M+ Me+...+ Mp 


c, L(x) (1.04) 


ty—™m » or (— X,—™M,, | AY A ok 
; Tho..-Kyt 
(m,—1;)!...(m,—rT,)! 


(1.05) 


+1 ms (x, + a Dine (a+1; l: 


rs "  m,! 

a +1, a+1,...,a,+1, +1; ky, kp,...,h,,1), (1.06) 
M(P+9) Ra 
I(p) 


Lauricella’s hypergeometric function of n+-1 variables. 


—™M,, —Mog,.... —My, —83 


where (p), = > —1 for all the parameters a, F, being 


After these preliminary remarks the following proposition can be 
proved: 


If we consider the general ret 


Sa) LP (x) = "Ye, Lee Pa -> b=, (1.07) 


whe re 


C, = ¢,(m, n, x, B), b, = b,(m, n, «, B), 


R(«) > —1, R(p) > —1, R(a+p) > —1, 
these coefficients satisfy the relation 
b,(m, n, x, B) = (—)"*"e,(m,n,B—m-+n,a+m—n). (1.08) 


[ prove this proposition by means of the Fourier method of deter- 
mining the coefficients b, and c, in the developments (1.07). It follows 
from the orthogonality of the Laguerre polynomials that 

= * 


[ e-eaa r(x) LP (x) Le P(x) de 


gs! 


P(st+ta+p+ t) 


0 


fr . »/ L(x) L(x) 


ds le Ly 8+ +B] dh 
D(s+ a + B+ 2 


dx’ 


(s = 0, 1, 2,..., m+). 
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Using the first representation (1.02) of L{*+(«) and integrating by 
parts s times, we get 


fe LypS+O+ pe [LP @)LP(a)] 2 
Ties de® 


Now, if we take into account the formula 


qT" LO (x : 
Fie) — (_y Leta), 


we obtain for the coefficients of (1.07) the expression 


l s ¥ oO 
—= —r 8+a+B (X+1) (> (B+s-—-r) 
C, Te+atP+l) >, (") | ety Diet? (a) Leto (a) daz. 
7 . (1.09) 


The last integral is known by a result of A. Erdélyi,* and it is 
equal to 
ae | n+B Jr (s-+o+B-+1) 
m—r]}\n—s-+-1 
x F,(s+a+B+1; —m+r, —n+s—r;a+r+1,8+s—r+1; 1,1). 
Next, it is easy to verify (for instance, by using the well-known 
integral-representation of Appell’s hypergeometric function F, in 
which the variable elements are equal to unity) that 
T'(c—b+b’) l(c’ +b—b’) F(a; b, b’;c,c’; 1, 1)/T(e) T(e’) 
— _ pie ee ee ;c’+b—b’,c—b+6’;1,1). (1.10) 
Putting a = c+c’—1 in (1.10), we deduce 
(—)P"Pe)P(c’) 
T(c—b-+-b’)l(c’ +b—b’)’ 


and by means of this relation we obtain the interesting formula 


F,(c+ce’—1; 6, 6’;c,c’;1,1) = (1.11) 


io 8) 


Peay ty L6)(a2) L(x) dx = (—)rta(PTH\(IT?) (RG, v) > —1). 
a q \Pp 


0 


* A. Erdélyi, ‘ Uber einige bestimmte Integrale, in denen die Whittakerschen 
M;..,-Funktionen auftreten’: Math. Zeits. 40 (1936), 693-702. See also K. 
Mayr, ‘Integraleigenschaften der Hermiteschen und Laguerreschen Poly- 
nome’: ibid. 39 (1935), 597. 
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Thus, by (1.09), 


C, = ¢,(m,n, a, B) = (—)m+"+8 > a 
r=0 


Afterwards, if we write the second expansion on the right of (1.07), 
we obtain at once 


b, = b,(m,n, a, B) = (—)* S > one Faia 


r=0 
= (—)™+"e,(m, n,B—m-+n, atm—n), 


which is the required relation between the coefficients of the expan- 
sions (1.07). Hence 
m+n 
L(x) L P(x) = > ¢,(m,n, «x, B)L+P(x) 
s=0 


m+n 


(—)mtn > c,(m, n, B—m-+n, a-+-m—n) =. (1.13) 
s=0 


If we put m = n in (1.13), it reduces to 


2m x 


2m 
L(x) LBMx) = Fc, L*Px) = > e,5 
s=0 


s! 
(R(a), R(B), R(a+B) >—1), (1.14) 


so that, if we expand the product of two Laguerre polynomials of 
degree m and of different parameters « and f in a series of Laguerre 
polynomials with parameter «+f or in a series in 2*/s!, the coeffi- 
cients in the two expansions are identical.* 

We can further remark that term-by-term computation of the 
Hankel transform of the two sides of the expansions (1.13) yields the 
integral equation 
[ erztP, 2 /(ey)} LI ce) LAP (a) dh 
0 

iis (—)mtme-vya+B) [(B—m+m) (y) [ia+m—n)(y) (R(a+8) > —1) 
recently proved by A. Erdélyi.t 


* The particular case « = 8 of (1.14) has been given by W. T. Howell, 
Phil. Mag. (7), 24 (1937), 1082-93. 

+ A. Erdélyi, ‘The Hankel transform of a product of Whittaker’s func- 
tions’: J. of London Math. Soc. 13 (1938), 146—54. 
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2. I established in the note quoted above* the expansion of the 
product of two Hermite polynomials of different degrees 


min(m,n) m\ (n 
. ~ wfeé € 
H,,(x)H,,(x) = > r12 \ Jpn) (2.01) 

r}\r 

r=0 
of which, by means of the Fourier method, I obtained a further proof 
suitable for generalizations. This method can be illustrated here 
by the case of the product of three Hermite polynomials of degrees 
m>n>p. Writing 
m+n 


H,,(2)H,(«)H, (2) =" SAM"? (x), (2.02) 


m 
s=0 
and remembering the orthogonality of H,(x) and the definition 
8(p—x* 
(er ) 
dacs 
we get for the coefficients of (2.02) the expression 
s! 9s Nar Alm mp) 
ail d| H,,(x)H, (x) H,(x) | 


dx’ 


= | A,,(x)H,(x)H,(x)H,(x)e™ dx dx. 


m 


© 


Now, 
ds(H,, H,,H,,) _ 3198 (" 


dx’ 


n\{p | 
NG) fe) a MH, ~§(2 )H,_,(a J 


a 


so that 


oO 
m\ {[n\ [p Pe 
Nar AQMD) 9 > ("\(\(t) [ H,(x)H,()H,(x)e-™ dev. 
i+j+k=er ide 
tepeating the same process for the evaluation of this integral, we 
get finally 


A (m,n,p) 
“*M+h+ p— or 


m!n! p! 


t)\9/\k]/\r—9]\r—2 
: : a (2.03) 
, 
k 
* E. Feldheim, loc. cit., p. 1. See also G. N. Watson, loc. cit. Formula 


(2.01), as I learn from Mr. A. Erdélyi, had already been proved by Dhar, 
Bull. Calcutta Math. Soc. 26 (1934), 57-64, equation (7). 
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If we put p = 0, (2.02) reduces to (2.01), because 1+-j = 2r, and we 
must therefore have «1 = 7 = r in wget Thus 


n 
Ane), = rial |"), 
r r 


which is identical with (2.01). 
A shorter proof of (2.02) by means of (2.01), communicated to the 
author by the referee of the paper, is 


H,, (a) H, (2) H, (2) 


= “— r! 2 (7) (anal )H,,(x) 
on > : r!p! aro(™) (") eo) (2) nn «ane 
ab ad p p 


a =F Ani ‘p— Seen (2), 


where 


4 oa 2S m!n! (m-+-n—2r)! p! eae 

Patio : sr! (m—r)! (n—r)! (s—r)! (m+n—r—s)! (p—s+r)!" 
This is in appearance slightly different from (2.03), but is really 
equivalent to it and simpler. On putting p = 0, s = r, we recover 
(2.01). 

3. In the second part of the present note I shall investigate the 
problem of integral representations of products of Hermite poly- 
nomials. For this purpose we make use of the Gauss transform of 
Hermite polynomials, expressed by the general formula* 


] r (y—x)*/a _@ am x a 
Tera) | e~ H, (3) dy =(1 -¥) Fal Toa}p (3.01) 


= 


« being a parameter given over the range (—0°0, +00), and R(a) > 0 
Two important limiting cases of (3.01) are 


oO 


x [ e-u-2 FT (y) dy = (2a)" (A=1;a—>1) (3.02) 


V7 


. 
—-x 


e | 
and — | e-u-iz!(2y)™ dy = 1H, (x) (A>0;a=1). (3.03) 

Nir 

— 2 
* W. Myller-Lebedeff, Math. Ann. 64 (1907), 388-416. See e.g. F. Tricomi, 

Ann. del Inst. H. Poincaré, 8 (1938), 111-49, where also further bibliographical 
indications are given, or E. Feldheim, Giornale Ist. Italiano Attuari, 8 (1937), 
303-27, and my previous papers quoted in this note. 
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Now, if we apply the Gauss transform (3.01) with A=1 to the 
identity (2.01), we get 


roa | e-u—-a¥la FT (x) H,, (x) dx 


min(m,n) 
"See oom tala} 8 
r=0 


On the other hand, we can establish other forms of the Gauss trans- 
form of products of Hermite polynomials as particular cases of a 
more general formula. The method applied here avoids the use of 
the expansion (2.01), and is based on the well-known generating: 
function 0 
ea? 
n 

If we employ the relation (3.05) in the variables x+y and x+2z, 
multiply their product by e-*/¢ (R(a) > 0), and integrate with respect 
to x over (—oo, +00) (the uniform convergence of the double series 
permitting the term-by-term integration), we get 


. | ela (a+y)H,(x+2z) dx 


= e%y+2Tz-+7") | eT la+2U+T)x dx 
= 
PoP Qo pptrTatr( |] — aq) +9(2a)r 
= (az) > = Sa = = H,| ] ai ela 
p=0 g=0 r=0 dilaiali v( ) v( 

The integral on the left of (3.06) is equal to the coefficient of 
tT” /m!n! in the triple series on the right. To select the coefficient 
of ¢“7™ in the last multiple series, we have to take p = m—r and 
gq = n—r, where r assumes all the values 0, 1, 2,..., min(m,n). We 
consequently have the required integral transformation 


f ela (x+-y)H, (a+z) dx 
(az) 


—@ 


min(m,n) n y z 
= (1—a)im+) > - r? )( |r| 3 =r a 
r (1—a) y(1—a) 


(3.07) 
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Putting z = y in (3.07), we obtain from this relation another 
expression for the Gauss transform of H,,(x)H,(x) already given by 
(3.04), namely, 


ale vila FT (x)H,,(a) dx 


(1—a)Hon+m) i (=) %0()(-) nel aa) | aaa 


r=0 
(3.08) 


In the limiting case a = 1, the two expressions (3.04) and (3.08) 
reduce to the same formula; if a #1, we deduce from (3.04) and 
(3.08) the identity 


ns (2 Vn) ("\ ans) 2H «(YH «(y)] = 0 
> l—a . r r m+n—2r\¥ m-r n—rb¥ ia 


r=0 


for all positive integral values of m, n, and for R(a) > 0. 
When we take a = 1 in (3.07), we have 


= min(m “nm n (2 )"-*(2z)"*— 
‘ (ety), n(@+z) dx a +o ” ") “y . 


We may suppose without loss of generality that m >n. Thus 


n 
m\{[n 
r1/ \ Jeym—r2ey 
rjy\r 
r=0 
n 
m \ (2yz = 
— yn! Qm yn +3 (2 - y — 9! gm ym—n Tm m)(—2yz), 
. S—f, Tr: 
r= 


and we obtain the following integral equation for the product of two 
Hermite polynomials 


] 
Nir 


- (3.09) 


eH, (a-++y)H,(a—2) dx = n!2mym—"L&"-™)(2yz) (m > n). 


m 


This equation contains many interesting special cases. For in- 
stance, if we put y = z = 0, (3.09) yields the orthogonality relation 
for Hermite polynomials. When we replace y and —z by ty in (3.09), 
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we recover the limiting case of (3.04) and (3.08) for a = 1 mentioned 
above 
co 
ae | er . 2 
— ea-iw HT (2), (x) dx = n! 2™(iy)™-" L'°"—-™(2Qy2)  (m > n). 


Vi 


~— (3.10) 
This formula was. first given by Howell.* Certain important par- 
ticular cases are worth mentioning here. 

The case m = n yields the interesting formula 


] 


N7r 


[ e—(2—-tuY’ H? (x) dx = n! 2"L, (2y?), 
<a 
which can be written also in the form 
fe 8) 


LA | H?(a)e-™ cos 2ay dx = n! 2"e-"" L, (2y?), (3.12) 


V7T 
ry 
00 


expressing the Fourier transform of H?(x)e-* by means of the 


Laguerre polynomials of parameter zero. Likewise 


] 
Nir 


[A n(&)H, (x)e-™ cos 2ay dx = (—)?n! 2” e-U’'y?P L2P)(2y?) 
~ (3.13) 
for m—n even (= 2p), and 


ao 


3 


Vir 


— (3.14) 
for m—n odd (= 2p+1). 
It is easy to give other expressions for the integral transform con- 


H,,(x)H,,(x)e-* sin 2ay dx —)Pn! 2m e-v'y2p+1 [2p + 1)( Dy?) 


mm 


sidered above taking into account the connexion existing between 
Hermite’s and Laguerre’s polynomials, namely, 
= 192 (— 4) (2 ae { D2m+1» (4) (4-2 
H,,,(2) = (—)™m! 22™ [0 (22), Hom4(%) = (—)™m! 22+ Dh?) (x*). 
(3.15) 


Starting with the generating function of Laguerre polynomials 


“ 


(1—-Rexp( — 3) = > ™E(x) (valid for |t] <1 
m=0 


we find quite simply the formula 


L(+ B(2x) = F L- V(x) LB-V(x) (R(a), R(B) > 
r=0 


m 


* W. T. Howell, loc. cit. 
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and, more generally, 


Lie+P)(a+-y) = pd Lie D(x )L‘B-Y (y). 


Taking « = 8 = 0 in (3.16) and using (3.15), we obtain 

> (7) Ho rtan-vly) = (—Yrnt 2, 2y), (8.18) 
r=0 

which, in virtue of (3.11), yields 


n 


] , tet = (7 
p—(x—tuyY FT? (x) da = (—4)" H.,,,.(y) Hg, — . (3.19 
[ ¢-@-wrHa(x) de = (—3) > () Y)Han-aoly)- (3-19) 


Vi 
ny 


Formulae obtainable in a similar manner from the preceding results 
are left to the reader. 

4. We conclude by discussing the inversion of the functional equa- 
tion, an integral representation for the product H,,(2)H,(x), which 
may be established by means of the well-known reciprocity formula 
of Fourier. We can deduce this inversion also from (2.01), applying 
the integral representation (3.03) of H,,(x). In fact, this relation 
transforms (2.01) into 


; © 
n . 9 
= [m\[n\ 1 I | \ ai ain 
1, (x) H(z) = > r! 2? — [ e—(y—tz) = dy (m>n). 
‘ r}\r} wr . a | 
“ (4.01) 
If we change the order of integration and summation, it is easy to 
see by the argument used in deriving (3.09) that the inversion of 
(3.10) is represented by 


r 


m 


H. (a A, (x) = <2 am fe (y+ ix} (iy) ym- n[(m—n)9 y? 2) dy (m > n). 


< (4.02) 
The formulae corresponding to (3.13) and (3.14) are thus of the 


same form 


1, (a) H,, (x)e-™ ait 2m — fe y* 


{cos ¢ 2a y\ ym—” Lm) (Dy?) dy 


(4.03) 


Nit \sin 2xy J 
according as m—n = 2p or m—n = 2p+1 (m > n). 
In particular, when m = n, we obtain the inversion of (3.12) 
@o 
' 1 ' 
H2(x)e-™* = —n!2" | e-Ycos 2xy L,,(2y?) dy. (4.04) 


Nit J 
—@ 
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On the other hand, replacing H,,,,,_.,(2) in (3.04) by the expression 


(3.03), and carrying out the summation, we deduce the integral 
relation P 


Gz | eros gle) Ha) de 
Na 


a _ a Qm e+ tyL—a)( jg, m—n [)(m—)(Da72) dx 
(1—a) ' 
“ve (m>n;0<a<l). (4.05) 
Other formulae arise when, in (3.07), we replace one or both of 
the polynomials of Hermite by the highest power of their argument. 


We obtain, after performing the required calculations, the relations 


8) 


fe #'lal 2(u—ia) "A, (x-+-v) da 


-— © 


min(m,n) 
— gim(]—q)in = 3% \" 1 m)\ (7 * ioe 
“Fe lo all 


r=0 
(0<a<\}), 


4 (az) 


hence, for a = 1 and u = 2, 


3 hy e-(e—iv)" (2x) (x) dx 
Ni 


min(m,n) m\ [n 
— {mtn = (—2y( \ Ja vo)". (Uv). 


rj/\r 
r=0 


Afterwards, 


s [ e-@'[ 2(u—ix) |"[2(v—ia) |" dx 


NTT 
min(m,n) m\ [n 
= 277 (7) lt) 0) 
Tre 


— 2 


© 
7=Q 


and, form = n, u = 2, 


— { e-(r—tuy'(2a7)2" da = > 2)'r n(| ) H2_(u) = (—)"A,,(u). 


In a similar manner we can obtain various other integral formulae. 
The computation of these relations may be left to the reader. 
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Added 8 December, 1939. I have obtained the more general 


relation 


/ [ e-w'laH] 
: (az) J m 


oe H (P42 a 
a n be 


(oe—ayae—an} “(e)(e) (0a) [ray 
(4.06) 


min(m,n) 


) 


containing all the formulae of §4, and also (3.07). In fact, by 
specializing the parameters a, A, w, and the variables y, z, we can 
reduce (4.06) to the results mentioned above. 





ON THE PRODUCT OF TWO ASSOCIATED 
LEGENDRE FUNCTIONS 
By W. N. BAILEY (Manchester) 
[Received 12 August 1939] 
1. AsouTt twenty years ago Dougall* gave an integral formula equi- 
valent to 


1 


(1 er ae mk) Po m( He) Pr m() du 


-1 


2™(k- 2m)! T(k—p-+-m + >) P(k—q+ m-- )P(k—r+m +4) 
vn(k—p)! (k—q)! (k—r)!P'(m +4) (k-+-m-+4) 





> 


(1.1) 
where p, g, 7 are positive integers with an even sum 2k, and any one 
of them is less than the sum of the other two. Otherwise the integral 
vanishes. 

In this result Ferrers’ definition has been used. Dougall proved 
his result for m > —1, but for non-integral values of m his definition 
of the function P?'(u) does not agree with that used by Hobson.+ 

An immediate deduction from (1.1) is that, with Hobson’s defini- 


tion of P7'(u), and p, g, m positive integers, 


2 4 y) > 
(4? —1)-#™ PP (ue) Po'(H) 
min(p—m,q—m) | ¢ 9y)\! 
om iy A r m A Pp ee A; oe (p T q— 2m —2 ) < 


Am ___ Dy\t 
r=0 Atta m-r (p+q 2r)! 


ei cee FR 
2p-+ 2q -2m—2r-+1 p+q-m—2r 


(3), j _(}—™m), 
es : ' , 


~(s+m)!’ r! 


When m 0, this formula reduces to the well-known Neumann- 


where AS 


Adams formula which expresses the product of two Legendre poly- 
nomials as a sum of such polynomials. When m-— oo, the formula 
gives the corresponding result for Hermite polynomials,t namely, 


a 2’ Anan oy (x) 
=, r!(m—r)! (n—r)! 
* J. Dougall, Proc. Edinburgh Math. Soc. 37 (1919), 33-47, formula (3). 
+ See $4 below. 
t W. N. Bailey, ‘On Hermite polynomials and associated Legendre func- 
tions’: J. of London Math. Soc. 14 (1939), 281-6. 


H,,,(«)H,, (x) = m!n! 


am 


(1.3) 
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due to Feldheim.* For general integral values of m, (1.2) expresses 


the product dP, (u) dP (u) 
du™ du™ 
as a sum of the type 
in( m,q—m) 
= — ss a OP 1 -m—2r(H) 
du™ 


, 
— 
r=0 
a formula which is useful in some problems of quantum mechanics.t 
So far as I know, no proof of (1.2) has been published. Dougall 
stated (1.1) without proof, and remarked that his method, though 
complete, was long and far from neat. Some years ago I used trans- 
formations of hypergeometric series to prove the formulat 


(p+m)! (q+m)! 
l im pm pm 7 
) p (u) q (14) 2™(p—m)! (q—m)! . 


9 
(* 
i 


q~-m A 


rm Ap _p An, (2p-+2q+ 2m—4r+-1\ 1, 
x 2? 4 m : ( : | a ) Pi gem-vt (1.4) 


9 9 9 or 
home 4 p+a+m-r 2p+2q+ 2m—2r-+1 


where again p, q, m are positive integers, and p > q+2m, q >m. 


When q—r is negative, (}),_, must be replaced by (—1)’-*/(),_,. 
This formula, like (1.2), reduces to the Neumann-Adams formula 
when m = 0. 


In the same paper I also showed that, when p, q¢, m are unrestricted, 


(? 4 Qu) Or me 1) . Qm—lemmin/rz yx 


4 (+m), P(p+m+r+10(q+m+r+1)0(p+qt+2r42) 
_ r! D(ptr+3)0(q+r+3)0(p+¢q+2m-+ 2r+2) 


r 


“~ 


b 1q+-m-+r+3 
, Teta tm ttt 8) (994 29-4 2m-+-4r+ 3) QM, smears): 


M(p+q+r+2) , (1.5) 
v0 


At the time when I gave these results, I failed to see how my 
method could be used to give (1.2). I now find that the same method 
can be employed to prove (1.2) and also to obtain some other 
formulae for products of associated Legendre functions. 

E. Feldheim, J. of London Math. Soc. 13 (1938), 22-9, (1.4). 
+ W.N. Bailey, loc. cit. 
- W.N. Bailey, Proc. Cambridge Phil. Soc. 29 (1933), 173-7 
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2. It is evident that, when |x| <1, 


(a, B; v3 2) F(a’, BY; y’; x) 
= ate" nlP In on mae = I- 8s e. | 
- y, 1—n--a’, 1—n—B 

We assume the relation 


yty’ = a+B+o'+P’, (2.1) 


and then the series ,F, on the right is Saalschiitzian, and can be 
transformed into a well-poised ,/, by Whipple’s formula* 





~[% Y, 2 —23 U—Y),(U—Z), 
al y |={ -y . 


u, Vv, W (u),(u—y—2)» 
x .F| ™ 1+}a, w—%, v—-%, Y, 2% —MN; 
la, v, w, 1t+z—n—u, 1+y—n—u, 1+y+2—u]’ 
where a= y+z2—n—u = v+w—a—l1. 
Taking x=1—n-y’, y=a, z=f, wu=1—n-—oa’, v= 
Y> ? Y> 


w = 1—n—Ff’, and expanding the ,f;, we get 


Fla, Biyit)Flel, Biya) = > > (— Eder ee ee 
(«-+B-+0’),-s(0-+B-+a’—1+ 2r)(y’—B')o(a),(B)(y ty +2—N), 
r! (n—r)! (y),(1—n—B’),(a’ +B),(a+a'),(a+B+a'+n), 


x” 4 


Writing n = r-+s, and summing with respect to s, we find that, if 
(2.1) is satisfied, 
F(a, B3 y; 2) F(a’, B’; y’; x) 
ii = ag B’)(«)-(B)-(y t+’ — 1) a, 2” 
"iar B-+a'-+r—Hly-+y’—H), 
<a fe be se te 
af'3 . 
y’ +r, atB+a’+2r, y+y’—1+r 


* See W. N. Bailey, Generalized Hypergeometric Series (Cambridge tract, 
1935), § 7.1 (1). 
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33 
The ,Ff; reduces to an ordinary hypergeometric function when 
y’ = a'+B, y+y’—1= a+a’, and then we obtain the expansion 


el % Bs . Tafa’, 1-8; 
a Peer | a’ +B | 
Po S (a),(B),(a’-+2B—1),(ata’+r), 
r «7! (a’+-B),(1+-a—B),(a+B+a'+-r—1), 


is i ia .— x]. (2.2) 


a+fB+«a’+2r 


In (2.2) all the series are well-poised, and can therefore be expressed 
in terms of associated Legendre functions. 


3. Now, if m and n are positive integers and m < n, it is known 
that* 
2m—2n(2n)! 7 ; ‘ 
P™(u) z"—™(y2— 1)" F(4-+-m, m—n; $—n; 2z-*) 
n!(n—m)! 


2 —m—2n(2n)! 


on! ——— ' —m—n;}—n;2~*), 


where z = p+./(u2—1). Using (2.2) with «= m—p, B = $+m, 
x’ = —m—q, we obtain (1.2). Taking «= —m—gq, B = }—m, 
m—p in (2.2), we obtain (1.4). This is equivalent to the method 
which I previously gave for (1.4). 


, 
x 


Again, when p, qg, m are unrestricted, we havet 


., T(n+m+l1), , 
¢ mw) — Qmemtin/y Ten ) (y2—1)¥mz—n-m-1 5 


x F($-+m, n+m-+1;n+8; 27") 


" 
. q-mgmetale | (ott 2), 2 ])-imz-n+m-1 x 
(n+) 


x F(4—m, n—m+1; +3; 27), 
where z is as before, and |z| > 1. 
Taking 


x = p+m+1, B = $+m, a’ = g—m-+1 in (2.2), we obtain 
Cf. E. W. Hobson, The Theory of Spherical and Ellipsoidal Harmonics 
(Cambridge, 1931), 114, (51) and (52). 

+ Hobson, loc. cit., 234 (67) and 196 (21). 


2695.11 


D 
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(1.5). Also, taking «a = g—m+1, B = }—m, a’ = p+m-+1, we ob- 
tain the formula 
(ptm + UP (q+m-+1) 
2__] smgym m = 2-m- 1 mrin 
(ut —1 Ou) OH) = Bo terete OT gamed) 
‘ oe 1m) P(p—m+r+1)0(q—m+r+1 \'(p+q-— ~m+r+4) 
r!D(ptr+ 30qt+r+ 30 (p+q+r+2) 
x (2p-+-2q—2m +4r-+3) pa q—m+2r+1(): (3.1) 
This result can also be obtained directly from (1.5) by changing 
the sign of m, and using the well-known formula 
emi n(n Be enrtg--™\) 
ee ~ T(n—m+1) 
4. As I remarked in my former paper, we can use the formula 
e-mri ein 


Pit() = ———_{ Qn'(u)sin(n-+-m)x— Q™ ,4()sin(n—m)z} 
(4.1) 


to obtain the product 

(uw? —1)-2" PP) Q7'(H) 
in terms of two series of functions of the second kind, in one of which 
the general term involves Q?. 94.n+2r+1(#), and in the other the general 
term involves Q”"" .m+2-(#)- A similar method can now be used to 
obtain the product (u2— 1)" PQ 


in terms of two series of functions of the sec a” kind. 

Dougall’s definition of P”'(~), when m is not integral, agrees with 
Hobson’s definition of P;™(u), except for a factor e””. Thus 
Dougall’s integral formula suggests finding a series for 


(u?—1)-"P> hes \P hea PM ), 


q+m 


where p and q are positive integers and m is not a positive integer. 
In this case (4.1) reduces to 


Demmi oj 
ae SIN M77 
2I—m ake —m 
I pt m( be ) ae — =» m-1(F)> 


7 
and so, from (3.1), we get 


9m»! a! 
2_] im pP m Pom —_ =" P-4: 
(HET pm ENE gtm(Y) = rp Om 1)P (q+ 2m-+1)P(4-+m) * 
in(p,q) . F ¢ 
, mR ff BL -+.m Bq-r.m(2p+2q+ 2m—4r+1 —m aa 
; rol Do e+n- r 2p+ 2q-+2m —2r+1 hil alias 
Vil i I we 
where B, m = I (3+m+r) Br = — I (3 T 8) . 
; [T(r+1) C(m+s+1) 
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This generalizes (1.2), to which it reduces when m is a positive 
integer. 


5. We now turn to the formulae* 


l p—l —km 2 —p p—l 
Po (4) = = (: ) ( r(— 2—m, —p; 1—m;— a 
em P—m)\e+1) w+ , ’ p+ 


i ag OE 1,p+1;1— 
- rtm) “Galen 


valid for R(u) > 0, and for general values of p and m. Using (2.2 
with a = —p—m, B = —p, a’ = p—n-+l, we get 
(u?—1)- Wl cers pet) = 
_- =. (—p),(—p—m),(—p—n),(1—m—n-+r),T(—p—m—n+r) - 
_ r! \(1—m+r)P(1—n-+r) 
xX (p +m+n—2 r) Pp m+n— *r(u), (5.1) 


and, by taking a = p—n+1, B = p+1, a’ = —p—m, we find an 
equivalent formula, remembering that 
P P(e) = : es (4). 
The formula (5.1) can also be derived from (1.5) by using Whipple’s 


transformation of Q7(cosh «) into P—%—4 (coth «). 


* Hobson, loc. cit., 210 (41). 








THE GENERALIZED BESSEL FUNCTION OF ORDER 
GREATER THAN ONE 


By E. M. WRIGHT (Aberdeen) 
[Received 11 May 1939] 


1. In a former paper (1) I used the method of steepest descents to 
determine the asymptotic expansion of the integral function 


x 
ait 


we 2. P(n+1)P(pn+B) 


for large z when p is real and positive. For p <1, ¢(z) has an 
exponentially small expansion in part of the plane; apart from this, 
the results for 4(z) for all p > 0 are particular cases of more general 
theorems proved by a different method in a later paper (2).* The 
zeros of ¢(z) lie near the negative half of the real axis. 

When —1 < p< 0, ¢4(z) is still an integral function, but it no 
longer satisfies the conditions of the general theorems of (2). The 
asymptotic behaviour of ¢(z) now presents certain new features; ¢(z) 


is exponentially small in one sector, exponentially large in two neigh- 
bouring sectors and, if —1 < p < —4, ¢(z) has an algebraic expan- 
sion in a fourth sector. If —1 << p < —}, the zeros of ¢(z) lie near 


the two lines 


arg2 = +42n(3p+1), 


while, if —} < p < 0, the zeros lie near the positive half of the 
real axis. 

As the asymptotic behaviour of 4(z) when —1 < p < 0 cannot be 
deduced from the general theorems proved in (2), it seems worth 
while to work out the complete asymptotic expansion of 4(z) in this 
case. This is my object here. 

I find it necessary to use both of my former methods. For part 
of the plane I express ¢(z) as the sum of two functions satisfying 
the conditions of the general theorems of (2). For the rest of the 
plane this method is still applicable, but leads to an algebraic 
expansion with zero coefficients;+ in fact, the expansion here is 


* The notation of (2) differs from that used here and in (1) ; in particular, 
the function ¢(¢) of (2) has no connexion with the ¢(z) of the present paper. 
{ See the first paragraph of § 5. 








ON A GENERALIZED BESSEL FUNCTION 37 


exponentially small and to calculate it I use the method of steepest 
descents. rn 
We write o = —p; since —1< p< 0, we have 0<o0<1. We 
write y = —z, choose argz and argy to satisfy 
—a7 < argz <7, —a7<argy<iam 
and write Y = (1—o)(ory)'¥1-0, 
Y, = (1—o)(o%e7*) #a-0), Y, = (l—o)(o%ze-7*)¥A-), (1.1) 
Hence 
Y,=Y (-—7< argz < 0), Y%=Y (0< argz <7). 
We use ¢ to denote a positive number (to be thought of as small); 
L and M are any two positive integers; and K is a positive number 
(not always the same at each occurrence) independent of x, y, z, and 
t, but depending on «, o, 8, M, and L. The constant implied in the 
O-notation is of the type K. We write R(w) and 3(w) for the real 
and imaginary parts of the number w. 
By the well-known asymptotic expansion of the T’-function, 
M-1 (—1)"4,, 


T'(ot+1—8) fe 
27% 1—o)t-M DT (t+ 1) a P{(—oyt-+B+}-+m}* 


1 
ri(i- 25 1.2 
+(na=oerp tiem) (2.3) 


when argt, arg(ot), and arg(ot+1—f) all lie between —z and 7: 
This defines Ay, A,, Ag,..., Ayy_1. We write J(Y) for the ‘exponen- 
tial’ asymptotic expansion 


M— 
Yi-fe-¥ (SA, ¥-"+0(¥-™)| (1.3) 


m=0 
and J(z) for the ‘algebraic’ asymptotic expansion 
“3 AB—1-2) 
> ae a +- O(28-1-Die), (1.4) 
fa ol (l+1)P{1+(B—l—1)/o} 

In the neighbourhood of the negative half of the real axis in the 
z-plane, ¢(z) has the ‘exponentially small’ expansion J(Y); in this 
region ®(Y) > 0 and so e-¥ is small. This expansion is also valid 
in the two neighbouring regions in which R(Y) < 0, where J(Y) is 
‘exponentially large’. In fact, we shall prove 


TueoreM 1. Jf \argy| < min{$z(1—o), 7}—e, then 


$(z) = L(Y). 
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If 0 < o < }, so that 7 < 37(1—o), the only region not covered 
by Theorem 1 is that in which 7—e < |argy| < 7, that is the neigh- 
bourhood of the positive half of the real axis in the z-plane. Here 


the expansion is the sum of two exponentially large expansions. 


THEOREM 2. If 0 <0 < } and |argz| < m(1—o)—e, then 

$(z) = 1(¥,)+1(%). 

We observe that, for 0 < o < }, R(¥,) < 0 and R(¥,) < 0 when 
argz = 0, so that J(Y,) and J(Y,) are both exponentially large in the 
neighbourhood of the positive half of the real axis in the z-plane. 

For o = }, however, R(Y,) = 0 and R(Y,) = 0 when argz = 0. 
Hence /(Y,) and J(Y,) are neither exponentially large nor exponen- 
tially small but roughly comparable with an algebraic expansion. 
In fact, an algebraic expansion does appear in this case. 

THEOREM 3. If o = } and |jargz| < m(1—o)—e, then 

$(z) = 1(%)-+ 1%) +S). 

When } <o <1, there is a region of the plane in which the 
expansion is algebraic. 

THEOREM 4. If } <o0< 1 and |argz| < }n(30—1)—e, then 

$(z) = J(2). 

This theorem and Theorem | still leave untouched the neighbour- 
hood of the lines argz = +4n(30—1) when }<o <1. Here we 
have a ‘mixed’ expansion. 

THEOREM 5. If } <o <1, jargy| < =z, and 

larg z-+42(30—1)| < a(1—o)—e, 
then ¢(z) = 1(Y)+J(z). 


2. Proof of a lemma by ‘steepest descent’ 


ot—t?+1l—c. 


? 


We write f= 


o—l 
g(v) is the regular function defined near v = 0 by 
{ 2—o (2—o)(8—o) , 

no) = I, g(v) = {1-+-- y - ——_— 2 

g( ) ? g( ) \ i 3 3.4 = 
so that f(t) = —4ho(t—1)*{g(1—t)P; 
and a,, is the coefficient of v?” in the expansion of 

(~ 1)m2™—41'(m+ >) 9)—BS q(47)¥-2m-1 
zo™—t+B(] —q)!-P (1—v)-Fig(v)} 
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in ascending powers of v. At a later stage we shall prove that 
A,, = 4», but for the moment we write 


; _;Mz=1 
1,(Y) = Yi-8e-3 7 4. Y-"+0(¥-™)| 
m=0 
and cannot assume that this is identical with J(Y). 
In this section we prove the following 


Lemma. Jf jargy| < (7—e)(1—o), then 
$(z) = L(Y). 

We denote by C a contour in the u-plane starting from infinity 
in the part of the plane in which —37+}e < argu < —4n—fe, 
passing round the origin in the counter-clockwise direction, and going 
to infinity in the part of the plane in which }7+ }« < argu < $n—He. 
Chen 1 — 

T(—on+f) 2m 
Fi 


yor—Ber dy 


and ¢(z) = - uP exp(u+zu’) du, 
FA 
the change of order of integration and summation being readily 
justified. 
If we now make the transformation* 
u = (cy)a-o% — = Yt, 


—_ 


‘ ad 2 
we have ut+zu? = — y(t—<) = —Y{1+/(}, 
Co 


—e 


1-BY1-Be-¥ 
and d(z) = — 2 | t-Be-Y 10 dt, 
r 


where [ is the transform of C in the ¢-plane. 
The ‘steepest descent’ curves through ¢ = 1 in the ¢-piane are 


given by f(y} = 0 (2.1) 
and Ri f()} = 0. (2.2) 
The first of these has two branches, one the positive half of the real 
axis and the other a curve encircling the origin once, its two ex- 
tremities asymptotic to arg? = +7 respectively. The latter curve, 
described in a counter-clockwise direction, we call I. The curve (2.2) 


* This transforms the ‘saddle point’ into ¢ = 1. 
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has also two branches, each the reflection of the other in the real 
axis. One, which we call [,, starts from infinity in the upper half- 
plane, encircles the origin in a counter-clockwise direction, and passes 
to infinity in the upper half-plane; the curve has no asymptotes, 
but on its extremities argt—» — $m and argt—> }z respectively as 
|t] > 0o.* [ is the other branch of the curve (2.2), also described 
in a counter-clockwise direction. The transforms of f, [, [, in the 
u-plane we denote by C,, C,, C3. 

Since |argy| < (7—e)(l1—o), we have |jargY| < z—e. We con- 
sider three cases, viz. 


—tr-+e <argY < }n-—e, (2.3) 
e<argY < m—e, (2.4) 
—mt+e <argY < —e. (2.5) 


When e« < jz, the cases overlap. 

When (2.3) is satisfied, it is easily verified that C, has the pro- 
perties laid down in the definition of C; hence [, may be taken for 
[. The transformation 


w = —i(to)'(t—1)g(1—t), ie. w? = f(i), 
transforms [, into the complete real axis in the w-plane, described 
in the positive direction. For sufficiently large MV and all real w, we 
can show that 


1—fy—B 2M—1 
¢ t dt = ¥ b,w"+O(w%), 
2ri(l—o)'-Bdw Hy" , 
a 
where 6,,,, = i py apart from trivial modifications the proof is 
m--s 


identical with that of Lemmas 5 and 6 of (1), and is of a familiar 
type. It follows from the above that 


x 


Q 


(M-—1 


o(z) = Y1-Pe-¥ | tee Tr “HH iJ wme-V dw +O | w2Me-ROW* day 
(m 


m=0 Py 


=e Yi-Pe-¥ (Sa Y-"+O(/Y m)\ 
\ ” J 
= L(Y), 
since R(Y) > KY}. 


* When o = },T, is the parabola 2y+1 = a. 
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When (2.4) is satisfied, C, may be taken for C and I, for fT. The 
transformation 


w = —el7(ko)'(t—1)g(1—t), ie. w? = if(t) 


transforms [, into the complete real axis in the w-plane, described 
in the positive direction. For sufficiently large M and all real w, 


ol-Pt-B dit _ sM>t ; 
nil —o)I 3 a 3 b,, wre—izin +4 O(w™), 
= (1 n=0 


and so 


“a, * » tri(2m+1) 


7) — Y1-8, -Y 
$(z) = Y1-Re o— 


fi. p2mei¥ w* dw + 


2) 


+0 f f w2Me-XY jw? au) 


(Mal 
= Yi-fe-¥| ¥ a,, Y-"+0(Y-*)} 
m=0 
= 1(Y), 
since 3(Y) > K|Y}. 
When (2.5) is satisfied, the proof is similar, C, and [, taking the 
place of C, and [. This completes the proof of the lemma. 


3. Proof of the theorems 
[t follows from the well-known result 


1 sin m(B— on)I'(on-+ +1 —f) 


T'(—on+f) — a 


that = | fenBid(ze-701) — ¢-nBid(zer01)}, (3.1) 
270 


-_ 


I'(on+ 1 —B) 1 rh 
rie+i) 


provided that ['(ot+-1—f8) is regular at each of the points 


t=0,1,2 


where 


For the moment we suppose this last condition to be satisfied; we 
shall show later (in § 4) that our results are still true when [’(ot+-1—) 
has a pole at one or more of these points. 
We take 
—n <argz <7, —n < arg(—2) <7 
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and write 
X = (1—«e)(o%x)"4-), ¢, = (—1T((+ 1— Blo} 


oI (/+-1) 
M 
1,(X) = 2nXt-Pex| > (—1)"4,,X-"4 O(X-¥)]}, 


m=0 


Jy(—2) = ¥ of —x) 8-1-4 O(aB-1-Dioy, (3.3) 


i=0 
The asymptotic expansion of d() is given by the theorems of my 
paper (2), since this function, unlike j(z), satisfies the conditions of 
these theorems. The « and x of (2) are replaced here by 1—o and 
o°(1—o)!-°x; X has the same meaning here and in (2). The formula 
(1.2) shows that condition A of (2) is satisfied in the appropriate 
region of the ¢-plane. We deduce from Theorems 4, 6, and 7 of that 
paper respectively that 
d(x) = I,(X) (3.4) 
when |argxz| < }n(1—o)—e, that 
d(x) = J,(—2z) (3.5) 
when |arg(—2)| < $a7(1+0)—e, and that 
d(x) = I,(X)+J,(—x) = 1,(X)+O(a26-2e) (3.6) 


when |arg2| < min{z, $7(1—o)—e} and |arg(—z)| < 7z. 

We now use these results in (3.1). We shall say that one asymp- 
totic expansion is dominant compared with another, if the second 
is of an order less than that of the error term of the first. Thus, 
by (1.3), when ®(—Y) > K|Y|, J(Y) is dominant compared with 
J(z) or O(28-Ye), and so 

I(Y)+J(z) = I(Y), 1(Y)+O(28-vle) = I(Y). 

First let us suppose that 0 < o < 4,80 that 7 < 3n(1—o). When 

largz| < $n(1-+-0)—e, 


1 
we have larg zet7°"| < $n(1+30)—e < m—e. 


Hence we can apply (3.6) to d(ze”™) and to d(ze-™”). In the first case 


arg 2z-+70 arg 2-+-7 


l—o l—o 


—n = argY,—7, 


arg X = 


by (1.1), and in the second 


,  argz—no argz—a i 
arg X = — = . = +a = argY,+7. 
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Hence, by (1.3), (3.1), (3.2), and (3.6), 


1 : oe . F 
¢(z) = 5; oP La(e"'Yq)—e- PL, (e-™"¥,)} + Of -?) 


= 1(%))+1(¥,)+ 0G?) 
For positive or zero argz within the interval considered, we have 
R(—Y,) > K|Y|, and so 1(Y,) is dominant compared with O(28-Hie); 
for negative argz, I(Y,) plays a similar role. Hence 
$(z) = 1(¥%,)+-1(%2), (3.7) 
4nr(1+c)—e. 
< n(1l—o)—e, we have Y = Y, and 


when |argz| < 
If e < argz 
, argz—7 argz+7 

cos arg Y,—cos arg Y, = cos 6? — cos 2824 

aang =— 


arg? sin” << Re 
l—o l—o 


= 2sin 
hence J(Y,) is dominant compared with /(Y,) and 
T(¥,)+1(%2) = 1%) = 1(Y). (3.8) 
Hence, by (3.7), 6(z) = 1(Y) when e < argz < }n(1+0c)—e. But, 
by the lemma of § 2, ¢(z) = 1,(Y) when 7—(z—e)(1—o) < argz < =. 
If « < 4n, the two latter regions overlap and so a,, = A,, and 1(Y) 
and J,(Y) are identical.* Hence ¢(z) = J(Y) when « < argz < 7, 
i.e. when —z7+e < argy <0. The same result can be proved 
similarly when 0 < argy < 7—e, and so we have Theorem | for 
Vv<co< 4. 
From this it follows by (3.8) that 
$(z) = 1(¥) = 1(%,) + 1%) 
when « < argz < a(1—o)—e and so (3.7) is true when 
0 < argz < m(1—o)—e. 
The proof for the corresponding interval of negative values of arg z 
is similar and Theorem 2 follows. 
Next let o = 4. The work proceeds as above except that we 
use the first part of (3.6) in (3.1) instead of the second. When 
in(1+o)—e, we have 


, - 1 , ; : . 
d(z) = 1(¥,)+-L( 2) + 5— {eid (—ze-7") —e-* Pt, (—ze™™)}. 
277 
* If for all |z| > K and for some fixed argz we have f(z) = I(Y) and 


f(z) = 1,(Y), then it follows that a,, = A,, and I(Y) and J,(Y) are identical. 
The proof of this will be obvious to the reader. 
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If « < argz < a(l—o)—e, then Y = Y, and J(Y,) is dominant com- 
pared with all the other terms. Hence, as before, [,(Y) = J(Y) and 


¢(z) = L(Y) (3.9) 


for —7+e < argy < 0. The proof of this result for 0 < argy < 7—e 





is similar and Theorem 1 follows. 
On the other hand, if |argz| < oz. then —ze+7% = ze*7(l-o% and 
arg |ze*7™-°%'| < w, Hence 


| fenbiy.( —ze-70) ey e~ Bi T / e ze7™7)\ 


5) 


210 
L-1 
aaa 1 pS ¢,28-A-Dle sin! np 4 mB—1—I(1—e)\ | gy 4p-1-1ye0) 
TT homed | Co } 
1=0 
= J(z) (3.10) 


by (1.4) and so ¢(z) = 1(Y,)+1(¥,)+J (2). 
By the use of (3.9) and considerations of dominance it is easily 
seen that this is true throughout the larger interval in which 





largz| < a(l—o)—e; 
this is Theorem 3. 
Finally let us suppose that } << o <1. If 





largz| < $n(30—1)—e, 
then arg(—ze+7")| = |argz+a(1—o)| < $n(1+0)—e 
and so, by (3.1) and (3.5), 
ble) = sferPidh (zer-on)—e-mPL (ze-m-o) 
= J (2) 
as in (3.10). This is Theorem 4. 
If 4n(30—1)+e€ < argz < }n(1+<)—e, (3.11) 
we have larg(ze-7%) | < 4n(1—o)—e 


and so —.erBid (ze-"2) — —)e*Bif,(eY,) ax J(¥,) = JF) 


a7 27 
by (3.4). Again arg(—ze7*') = arg(ze-74-9) and so 
—4n(3—50)+e < arg(—ze™") < 4n(30—1)—e. 
by (3.11). Hence 


jarg(—ze7")| < 34max(30—1, 3—5c)—e 


< }n(1+oc)—e 
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since }< ao <1. 


and so, by (3.1), 
$(z) = 1(Y)+O(e8-Ne) = I(¥), 
since J(Y) is easily shown to be dominant. The same result for 
—4n(l+o)+e < argz < —4n(30—1)—e 
is found similarly. Hence ¢(z) = J(Y) when 
4n(1—o)+e < jargy| < $a(1—o)—e. 
Combining this with the lemma of §2, as we did for the case 
0 < o < 4, we obtain Theorem 1. 
Without loss of generality we may suppose that 
e < 47min(30—1, 1—o). 
We can deduce (1.5) from Theorems 1 and 4 by considerations of 
dominance except when 
4nr(30—1)—e < |argz| < $n(30—1)+<. 
Let us suppose that 
0 < 4n(30—1)—e < argz < 4n(30—1)+€. 
It is easily verified that the conditions of (3.6) are satisfied for 
d(ze-7*) and those of (3.5) for d(ze™*). Hence 


em Bi 


= 1 ‘ ; : ' 
T,(e*'Y,) + —{e"P'(—2e-7) —e- "PJ, (—ze7")} 


$(z) = 


Qr1 
= 1(¥,)+I(2) = 1(¥)+J(e) 
as before. The proof for the corresponding negative values of arg z 


is similar. 


4. Removal of a restriction 

At the beginning of §3 we assumed that I'(ot+-1—f) is regular at 
each of the points t = 0, 1, 2,.... If this assumption is false, this 
function has poles at a finite number of these points. For simplicity 
let us suppose that there is just one such pole at ¢ = m,, so that 
on,+1—f = —l,, where J, is a non-negative integer. The argument 
can be immediately extended to the case of more than one such pole. 

If we replace the term in x” in d(x) by the residue of 

al (ot+1—Bf)(—2)! 
sin wt P(¢t+-1) 
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at its double pole at t = n,, we have the same asymptotic expansion 

for d(x) as before, except that we must put c, = 0.* This residue is 
b,x"+-b,x™ log(—z), 

where b, and 6, are independent of x and 


L Po. (— ] jf 

"2 GP(n, +P G+)) 

Hence the change in our results, if any, is always of ‘algebraic’ order 

and so may be ignored when our expansion is exponentially large. 

Since the work of §2 still holds good, there is no change in the 

exponentially small expansion ¢(z) = /(Y) valid in the sector en- 

closing the line argy = 0. Hence our results are unchanged except 

possibly when J(z) is of importance, that is, when } < o < 1 and 
jargz| < 4n(30—1)+€. 

Here we have 

gm 


1 . . ; ; 
— ___ SpmBid(opo—no0i)__ p—mBi (sp7ci = — 
$2) = 5 feMaleet0)—e-*dlcer™t} + 


—slb gita(em(B—ny 0) _ e-7(B—nyo)) 4 


+b, zmfeB—n Nog (ze71—2)) — e—wi(B-n Nog (ze-71—2))¥] 
If we suppose « < }n(1—a), we have |argz| < oz. Since we have 
to put c, = 0, the asymptotic expansion of the first term on the 
right is altered by the subtraction of the term 

2(B-1-l)lo 
oT (1, +1) {1+ (B—1,—1)/o} 

Hence the asymptotic expansion of 4(z) has the additional terms 

2M 2(B-1-h)lo 


P(my+1)P(—on, +f) of (,+1)P{1+(B—1—1,)/o} 


a —{(by+bs log 2)sin (B—n, o)+-m(1—o)b cos m(B—n, o)} 


o 


zm (12+ (—1,1—o)P (+ DP +D)} 


ees ml | a 
D(n,+1)0(¢,+1) 
= Q, 
and so our results are unchanged. 


* See the paper (2), especially the remark at the beginning of § 1.3 and the 
definition of H(—) in § 2.1. 
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A further result in a special case 

If jargy| < }n(1—o)—e, the condition of (3.5) is verified for both 
d(ze-7%') and d(ze™"), since —ze+™* = ye+™%, Hence, by (3.1), 
arg(—ze+7*")| < @ and so 


¢(z) = 5 ern (ye-72%) — e-*Bi J (ye7™™)} 


L-1 
=) > cylHe sin a(l+1)+ O(y8-Die 
"& 


os O(y8-1-Lie) 


for all L. In fact, in the lemma of §2 we succeeded in replacing this 
order result by an exponentially small expansion. 

Similarly, if the coefficients in J(z) vanish for all / greater than 
some I,, it is possible to replace the error term O(28-!-Die) by 
an exponentially small expansion. These coefficients vanish only 
if T'{1+(B—t—1)/o} has poles at t=1,+1, J)+2,..., that is, if 
| +(8—Il—1)/o is integral for alll > J). If this is so, 1/o and B/o are 
integers. If o < }, J(z) never appears (Theorems 1 and 2), while, 
if o = 4, J(z) only occurs added to an expansion whose error term 
is at least of algebraic order (Theorems 1 and 3). Hence we have 
only to consider the case o = }. 


THEOREM 6. If o = } and 28 is an integer, then the coefficients in 
J(z) all vanish after a certain 1 =1,. If H(z) is the finite sum con- 
sisting of the non-zero terms of J(z), but excluding the error term, then 
H(z) is a polynomial and 

$(z) = H(z)+-1(Y) 


when jargz| < ?7—e. 


Let oc = } and b = 28 be an integer. For all z 


A(z 
, P(n+1)0(—43n+}0) 


@ 
2m gam 


zr (2m+1)P(—m-+4b yt 2 (2m+-2)'{—m-+4(b—1)} 
= F,(z*)+2F,(z*), (5.3) 


where F,(x) and F,(x) are integral functions of x. 








48 ON A GENERALIZED BESSEL FUNCTION 


When jz7-+e < |argz| < ?7—e, I(Y) is dominant compared with 


= 


H(z). Hence (5.1) follows from Theorem 1. Next let 


jargz| < }a+e < fn. 


If we write Y’ = (1—o)(0%z) "4-0, 


we have Y’ . Y = 3, and Y’ ‘ex*mi, Hence 


_ M-1 
Y= Feet S A, ¥’-"+0(¥'-™)} 
m=0 
a et7t(1—b) J ( Y) (— 1°17 (Y). 

If we apply Theorem 1 to ¢(—z) we have ¢(—z) = I(Y’) since 

argz! < 4r. 
] 
If 6 is odd, = 0 
b is odk i—m+i—)) 
when m > 3(b—1) and so F;(z*) is a polynomial in z?. Hence, by (5.3), 


$(2)—d(—z) = 22 F(2*) = P(e), 


a polynomial in z, and 


= P(z)+I(Y). 
If b is even, F,(z?) is a polynomial, 
$(z)+¢(—z) = 2F,(z7) = P(z), 
a $(2) = Ple)—$(—2) = P@)—1(Y’) 
= P(z)+1(¥), 
which is (5.4). 
Comparing (5.4) when argz = 0 with the result of Theorem 5, viz. 
¢(z) = H(z)+ O(z8-1-De) 
for all L, we see that P(z) is identically equal to H(z). Hence 
Theorem 6 follows from (5.4). 
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SOME INVERSION FORMULAE 


By M. M. CRUM (Ozford) 
[Received 22 November 1939] 
1. A Fourrer kernel may be described as a function k(x) such that 


a formal solution of the integral equation 


oO 


f(x) = [ k(wy)g(y) dy | (1) 


0 


g(x) = [ ef) dy. 
0 


The theory is given by Titchmarsh;* it covers the examples given 
below. The theory is stated in terms of k,(x), the integral of k(x); 
the general formula for k,(2) is 


t+io 
[ aos oe? ds, (3) 


4-—io 
where &(s)b(1—s) = 1. (4) 
2. If, in (3), o(4-+it) = efasinht (a > 0), 
then, by Cauchy’s theorem, we can move the contour of integration 


a distance }z to the left. This gives 


8) 


x }7r+4 e—acosht— itlogz 
k, (2) | dt. 


int 


—- 2 


Differentiating we get 


k(a) ll e-4 cosh/—itlog x dt 


1 _ 
— — tt K tog 2(4)> 


K,(z) = 4 [ e-zcosht-vt d¢ (R(z) > 0). (5) 


~ 


- @ 


* E. C. Titchmarsh, Theory of Fourier Integrals (Oxford, 1937), Chap. 8. 
G. N. Watson, Theory of Bessel Functions (Cambridge, 1922), § 6.22 (7). 
3695.11 E 








50 M. M. CRUM 
If we put 
z2=é6, y=e, 2f(xz)= 


then (1) and (2) become 


F@ = J er€mcgnle) tn) an 


eit +m K. e+ pa )F(m) dy. 


If, as is usual in the theory referred to, f(x) is of L*(0,00), then F(&) 
is of L?(—00, 0). 


3. As an example of the transformation let 


eb™G(n) = Kyz,,(b)e-”® = (b > 0), 


so that — e-#7€F(¢) = | Ken (@)Kugag (be? dy. 


Now, if p(t) and q(t) are of ——— oo), and 


] 1 


P(n) = | p(tye™ dt, Qn) = — [ q(t)ei dt, 


4/(277) 


\ (277) 


ty 
- © 


oe ms 


then* [ P(n)Q(n)e-*” dn = [ p(t)q(@—t) dt. 


Let p(t) a e-acosht+igt q(t) ios e—bcosht+ift. 


then, by (5), 


Pon) = |(Z\Krel@ — Qn) = | (E) Kgs 


a) 
and so e-int F(é) a } [ e-voomhsy itg—bcosh(O—t)+iLO-4) dz. 


Let c = (a?+ 2ab cosh 6+5?)!, 
asinh @ = csinha, bsinh 6 = csinhB (a, B real). 


Then a+b cosh 6 = ccosh, atpB = 0. 


* Titchmarsh, ibid., Theorem 64. 
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Hence 


e-ivé F(£) —— i e—ccosht—B)+iE—Lt+il0 dt 


= Kyg_o(c)e®P +a, 
This holds by analytic continuation for | /(@)| < 7, since 
K;,(a) = O(e*7") 
and the integral defining F(£) is uniformly convergent for 
|1(0)| < 78. 
If we put 
6=in--C) (0<C<7), a = 1A, B = 7B, 


we have 


a 


1 [ K. (a) Ky¢4)(b)er—™ dy == Kye_p(c)e 8-44, 


i ) 
7. ilE+n 


=~ @& 


where A, B, C are the angles of the triangle whose sides are of 
lengths a, b, c. 


4. The formal solution of (1), when k,(x) is given by (3), but (4) is 
not satisfied, is - 
g(x) = | h(ay)fy) dy, (7) 
0 
x“ $+io 
x x*p(s 
where h,(z) = [aw dy = = | ———- da, 
0 


ery 


4-—io 


and (s)p(1—s) = 1. 
Taking o(4+it) = e-tacosht (a > 0), 


we have $(}-+it) = efaeosht 
Ps - a ‘ 
k(2) — — Kitog2(4), h(x) = Kjiog2(—*4). 


Making the substitutions (6), we obtain from (1) and (7) 


Kyg+p(04)G(n) dn, 


Ky¢+5)(—1t4)F(n) dy. 
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5. If in (3) we put 


. eiasint (a + 0), 


@ 


re eta sin/—itlog x dt 


we have ——— + 
(j—it)(3—it) 


- 
xv 
Da 
a7 
. 
0 


x 
We may put GaN = + 
n oc 
and integrate term by term; the calculus of residues then gives 
[log zx] 
k(x) = > J,(a)(xein—el), 
n x 


flog x] 
from which k(z)= > e”J,(a), 


n x 
for logx not an integer. This is a step-function (because y(4-+-7#) is 
g ; 2 
periodic) and not an integral; consequently (1) and (2) must be 
replaced by ; 
: rf «# 
fw) =— | g{%) dkyly), 


0 


ate) = = | s(2) ata. 


"oa 
0 


The substitutions (6) now give 


F(é)= > J,(a)G(n—6&), 


G(é) = s J, (a) F(n—€), 


which is easily verified for suitable F(é). 








ON THE FREE PATHS IN THE SUBSTRATUM 
AND THE GRAVITATIONAL THEORY OF THE 
ORIGIN OF COSMIC RAYS 


By G. J. WHITROW (Ozford) 
[Received 14 November 1939] 


1. Introduction and summary 

As a result of his recent cosmological investigations Milne has put 
forward the view that cosmic rays may have been produced by the 
acceleration of freely moving particles in the gravitational field of 
the universe.* 

For comparison with the universe of extra-galactic nebulae Milne 
began by constructing a uniformly expanding substratum or con- 
tinuous distribution of fundamental particles in uniform relative 
radial motion. According to the observer associated with each funda- 


mental particle the substratum is isotropic and obeys the cosmo- 


logical principle, in virtue of which each observer, using congruent 
measures and plotting his observations in Euclidean private space, 
describes the system in the same manner. It follows that, according 
to each observer, the substratum has a distribution of particle- 
density ndadydz, where n, the number of particles per unit volume 
near (2, y,2) at epoch ¢, is given by 
Bt 

c3(t2—P2/c2)?’ 
P being the vector (x, y,z), B a constant of zero physical dimensions, 
and c the speed of light. The velocity V of any particle of the system 
at position P at epoch ¢ is given by 

V = Pie. 

In the presence of this substratum the equation of motion of a free 

particle moving with velocity V through position P at epoch ¢ has 
been found, by purely deductive methods, to be of the form} 


dV Y Y 
= (P—Vi)+ G(é), (1) 


* Milne, Relativity, Gravitation, and World-Structure (1935), chap. 12; 
referred to here as WS. 
+ Milne, WS, § 101. 
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where 

X = #?—P?/c?, Y = 1—V?/c?, Z = t—P.V/c? 
and = FAY. 
This equation involves an arbitrary function G(é) and, from the 
manner of its derivation, is the equation of motion of a free particle 
moving in the gravitational field of the most general system of 
particles obeying the cosmological principle, according to the funda- 
mental observers of a uniformly expanding substratum. It has 
been shown by Milne that it is possible to supplement such a sub- 
stratum with a statistical system of free particles distributed accord- 


ing to the law N(P, V, t) dadydzdudvdw, 


where N(P, V,t) = i. 

so that the cosmological principle is still satisfied. Hence the acce- 
leration of a free particle moving in such a system must be included 
in formula (1). A system comprising a statistical system superposed 
on a substratum is called a ‘mixed’ system and a ‘pure’ substratum 
is thus the limiting case of a ‘mixed’ system when the statistical 
component vanishes. Consequently, to determine G(£) uniquely for 
a substratum we seek first its dependence on y(€) in the case of a 
‘mixed’ system and then the limiting form of G(é) for vanishing (&). 
In this way it has been found that for a ‘pure’ substratum G(é) = —1, 
and the equation of motion of a free particle is, accordingly, 


dV Y 
-—- —(P—Vi) =. (2) 

Milne has, however, considered in great detail the equation of 
motion of a free particle in the presence of the most general type of 
‘mixed’ system, where ¥(€) is an arbitrary function of € and G(&) is 
therefore also arbitrary but related to %(£) in a definite manner. He 
found that any free particle moving in such a system attained the 
velocity of light in a finite time. 

A ‘pure’ substratum can be regarded as a model of the universe 
in which all nebular condensations have been smoothed out. Its 
analysis can therefore be expected to yield insight into the general 
nature of world gravitation. A ‘mixed’ system, however, admits 
condensations but does not allow for any clustering effects of groups 
of nebulae or other irregularities of distribution. Consequently, in 
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analysing such a system in order to gain deeper insight into the nature 
of world gravitation we must proceed with extreme caution. How- 
ever, if we can legitimately regard a ‘mixed’ system as a model of 
the universe of extra-galactic nebulae, then, following Milne, we 
must conclude that a free particle at large in the universe attains 
the speed of light at a finite epoch. From this result Milne has 
derived a gravitational theory of the origin of cosmic rays which 
accounts in an entirely unforced manner for the following facts: 

(1) the corpuscular character of the primary radiation; 

(2) the high energy of the primary corpuscles; 

(3) the isotropy of the primary radiation and its independence of 

sidereal time. 
Unfortunately, as Blackett has pointed out, this theory is subject 
to a grave objection, for it would appear to involve equally the 
existence of uncharged as well as charged particles. So far these 
have not been found.* 

Since, from the experimental evidence, it is extremely unlikely 
that cosmic-ray phenomena can be explained by purely gravitational 
arguments, it appears that a ‘mixed’ system constructed in accor- 
dance with the cosmological principle does not yield verifiable results 
about world gravitation but is useful only in so far as it is required 
in order to obtain the properties of a ‘pure’ substratum as those of 
a limiting case. It remains, however, to prove that the same objec- 
tions do not arise in this limiting case. Many of the formulae asso- 
ciated with a ‘mixed’ system cease to be valid when G(é) = —1. It 
is therefore the object of the present note to re-examine the free 
paths in a ‘pure’ substratum and to prove that there are no pheno- 
mena associated with them which we are obliged to interpret as 
cosmic rays of a gravitational origin. Incidentally, a new least-action 
formula will be derived, since the general formula obtained by 
Walker,+ for the free particles in a ‘mixed’ system, breaks down 
when G(é) = —1. 


2. The free paths in a uniformly expanding substratum 
According to any fundamental observer O, keeping time ?¢, the 

substratum expands uniformly from O as point origin at zero epoch. 

At any subsequent epoch the system occupies a sphere of radius 


* Blackett, Halley Lecture (1936), 14. 
+ Walker, Proc. Royal Soc. A, 147 (1934), 484. 
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r = ct. Consequently, the path of a free particle which is not a 
photon, when extrapolated into the past, will be found to originate 
at O at ¢= 0 and at any subsequent epoch ¢ the particle will lie 
inside the sphere r = ct. Moreover, since it is not a photon, its speed 
will, in general, be less than c. 

At each instant after t = 0 a free particle passes a unique funda- 
mental observer. Relative to any one of these observers the motion 
of the free particle is rectilinear and is governed by the equation 

dV ,, 1—V2/c? 
—— = —(r— Vt) ____., (3) 
dt t2?—r?/c? 
where r is the radial coordinate and V is the radial velocity. If the 
clock t, kept by the observer, is regraduated to 7, where 
t = t, log(t/t,) +t, (4) 
then, with the usual conventions, the new radial coordinate A and 
epoch 7 assigned to an event at distance 7 and epoch ¢, according 


to the original clock, are given by* 


t 
A= 4etylog’ 


t- 
t = 3t, log 


dd Vt 


It follows that —= ——s 
dr t—Vri/c? 
For a free particle this r-scale velocity is constant.} For, if 
Vy = dafdr, 

 rt+het 
t+rV/c?’ 
dV (1—V3/e%(Vt—r) , df #—r2/c® | 


and so = Fe — ole 
dt (t-+-rh,/c2)? 


a . “ 1— V2/c? 

Since, from (7), : - = = 

(t+rK,/c?)? 

it follows by comparing (3) and (8) that dV,/dt = 0, since r? < c#?. 
Hence J, is constant and so 


then J 


sh - 0, (9) 
dr? 


* Milne, Proc. Royal Soc. A, 158 (1937), 344. 
tT Referred, of course, to an observer whom it passes at an epoch after ¢ = 0. 
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throughout the entire career of a free particle relative to a funda- 
mental observer whom it passes at a finite r-epoch.* 

From (7) we see that when r = 0 either ¢ = 0 or else V = Ky. Con- 
sequently, Vi, is the t-scale velocity of the free particle at the instant, 
other than ¢ = 0, when it passes the observer. Moreover, since 

r/t+VY r 

14+Vr[c%t~ t’ 

it follows that V > r/t and hence from (8) we find that dV /dt > 0. 
Thus V is a monotonic increasing function of ¢. When ¢— 0 and 
when ¢ + 0, then V ~ r/t, and in both cases} we find that V? > c*. 
Since V is a monotonic increasing function, we see that V > —c as 
i+ 0Oand V +cast—oo. Consequently, the career of a free particle 
according to a fundamental observer O keeping ¢-time, whom it 
passes at an epoch later than ¢ = 0, can be represented by the 
following diagram: 

>V=-c 


O 





t—0 |v =o. 





O 
V=J 
Fie. 1. 

The general features of the trajectory of a free particle relative 
to an observer whom it does not pass at any epoch subsequent to 
ti = 0 can be readily deduced. The path described is a plane curve 
emanating from the observer at ¢ = 0 and tending to infinity in a 
definite direction as t-> 00. The speed of the particle as t > 0, and as 
t + 00, is c, but at all other epochs the speed of the particle is less than c. 


3. The principle of least action 

The 7z-scale velocity of a free particle relative to any observer 
whom it passes at an epoch other than ¢ = 0 has been shown to be 
a constant, and the path is a geodesic in the observer’s private space. 
Hence, according to any such observer, the particle moves so that 


8 | dr = 0. (10) 


* Milne, Proc. Royal Soc. A, 159 (1937), 186, deduces this result from a 
least-action formula. 

+ This result has also been proved by T. Lewis, Phil. Mag. 20 (1935), 1101, 
and by Milne, Proc. Royal Soc. A, 154 (1936), 51. 
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But 7 is an invariant for all fundamental observers.* Consequently, 
relative to any fundamental observer the particle moves so that (10) 
is valid. Hence transforming to ¢-scale coordinates we can deduce 
a principle of least action. Since the r-epoch assigned to an event 


\v|=c : 
4 


Fig. 2. 


at position P at epoch ¢t, on the ¢-scale related to the 7r-scale by 


formula (4), is given by 


it follows that 
Writing, as usual, 
X = #—P?/c?, Y = 1—V2/c?, 
Z = t—P.V/c?, = XY 
Zdt Yidt 
—__ — 4 
oy = bet 


Now & is known to be constant throughout the career of a free 
particle.t Consequently from (10) it follows that 


we find that dr =t 


* Milne, Proc. Royal Soc. A, 159 (1937), 177. 

+ Milne, WS, § 160. This result may be verified in the following manner. 
From equation (7) above we find that €' = (1— V2/c*)-*. Hence is constant 
for a free particle relative to an observer whom it passes at an epoch after 
t= 0. But é is invariant for all fundamental observers. Consequently, é is 
constant along the path of a free particle relative to any fundamental observer. 
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ds 
8 | i (13) 
where ds = Ytdt = {dt?—(dx?+-dy?+dz?*)/c?}4. 
Hence, using‘t-time and Euclidean private space, the free paths are 
geodesics in public space-time of Riemannian metric 
ot = ee (14) 
P—(x*+y?+2*)/c? 
4. Clock regraduation and velocity invariance at the 
observer 
For comparison with observation we are primarily interested in 
the behaviour of a free particle as it passes an observer. We have 
already proved that the r-scale velocity of a free particle according 
to any observer whom it passes is constant. Moreover, this velocity 
V, is the same as the ¢-scale velocity of the free particle at the instant 
of passing. This result is a particular case of a general theorem: 
The velocity of any particle as it passes an observer is independent of 
the observer’s choice of time-scale. 
For suppose the observer O regraduates his clock from ¢ to 7’ where 
T = p(t) (15) 
is a single-valued, monotonic increasing function of ¢. Let a signal 
leave O at time t, by the t-clock, be instantaneously reflected at 
a neighbouring point, and return to O at-time ¢+-5¢. Let the corre- 
sponding epochs on the 7’-clock be 7’, 7+67'. Then, according to 
the ¢-scale, the radial distance of the neighbouring point is 
81 = 4ct. 
Similarly, according to the 7'-scale, the radial distance of the point is 
5L = 4c8T = fep'(t) de. 
Hence bL = y(t) al. (16) 
Now consider any particle passing O at time ¢. At the subsequent 
conventional ¢-scale epoch ¢+At let it be at distance Al from O. 
Then, according to the 7’-scale, at the subsequent conventional epoch 
T+AT the particle will be at distance AL, where 
AL = W'(t) Al 
AT = w(t) At J 
AL Al 
AT At’ 


(17) 


Hence (18) 
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from which result we deduce, in the limit, that the velocity of any 
particle, as it passes the observer, is independent of the particular 
time-scale chosen by the observer. 

The velocity of a free particle in the substratum is the same 
according to each fundamental observer whom it passes, if these 
observers adopt congruent 7-scales. Hence, whatever time-scales the 
observers choose, the velocity of a free particle in the substratum is 
the same for each observer whom it passes at the instant of passing. 
This is a result of great generality because it is true even when the 
observers concerned choose clocks which are not congruent. 

In the ¢t-scale the invariance of the velocity of a free particle 
relative to the observer whom it passes at any instant provides us 
with an immediate verification of the known result that, throughout 
the career of a free particle, €! is constant. Moreover, we can readily 
deduce that this expression is invariant under a new type of trans- 
formation. For, if the fundamental observers regraduate their clocks 


in congruent fashion by the relation 
T = kf, (19) 


they continue to describe the substratum as uniformly expanding 
from a point origin at zero epoch.* Hence, if each observer adopts 
Euclidean private space in association with the new scale of time, 
then in the new coordinates the expression €! is again invariant. But 
é} = (1—V?/c?)-!, where V is the velocity of the particle concerned 
(free or constrained), relative to the observer whom it is passing at 
a given instant. This velocity is the same for both scales 7' and t. 
Hence, the value of & which any fundamental observer, adopting 
Euclidean private space, assigns to any particle at a given epoch is 
independent of the observer’s particular choice of t-scale from the family 
given by (19). This result is of interest in view of the identification 
of m&! as the measure of the mass of a particle.t Incidentally, we 
note that it is not true that &! is invariant for a transformation of 
time-scale of the type (19) alone, because in general such a trans- 
formation will not preserve the Lorentz group. We must associate with 
each time-scale an independent choice of Euclidean private space. 
We conclude this section with one further application of the 
theorem that the velocity of any particle as it passes an observer 


* Milne and Whitrow, Zeits. fiir Astrophys. 15 (1938), 284. 
+ Milne, Proc. Royal Soc. A, 154 (1936), 31. 
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is independent of the observer’s choice of clock. In transforming 
from ¢ to t by (4) we have proved that the 7-scale radial velocity 
of a particle moving with radial velocity V at distance r at epoch t, 
according to the t-scale, is given by 

_ V—rft 
~ 1—Vr/et’ 
Milne has remarked that this formula is similar in form to the Ein- 
stein formula for the composition of the radial velocities V and r/t 
but, he adds, ‘it has a radically different origin’.* In point of fact, 
however, this formula can be immediately derived by means of 
Einstein’s formula. For consider a particle moving radially with 
respect to a fundamental observer O. At any instant it passes some 
fundamental observer O’, say. Let O, O’ adopt congruent 7-scales. 
Then the velocity of the particle with respect to O, when it passes 
O’, is the same as its velocity with respect to O’. But its velocity 
with respect to O’, at this instant, is independent of the choice of 
time-scale made by O’. If, according to O, the t-scale velocity of the 
particle at this instant ¢ is V and if O’ is at distance r then, since 
O’ has the t-scale velocity r/t, it follows by Einstein’s velocity- 
addition formula that the velocity of the particle at this instant, 


(20) 


v 


according to O’, is V—rit 

1—Vr/ct 
But we have already remarked that the velocity of the particle at 
this instant, according to O’, must be v. Hence (20) can be estab- 
lished directly by means of Einstein’s formula. Incidentally, it is 
worth pointing out that whereas the formulae 


t c 
= ct, log = v 
—r?/e2 


9 
9 


#2 
+ = tt, log 


+b 


involve ¢, explicitly, the formula 

_ V—rit 
aes 1—Vr/c*t 
is independent of t). This point is significant in view of Milne’s 
suggestion that at each epoch the appropriate 7-scale is given by 
choosing t) as the reading of the ¢-clock at that epoch. 


* Milne, Proc. Royal Soc. A, 158 (1937), 344. 
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5. Free particles and cosmic rays 

We have seen that, except as t>0 and as too, the #-scale 
velocity of a free particle in the substratum is always less than c, 
whereas in a ‘mixed’ system a free particle attains the velocity c at 
a finite epoch later thant = 0. Hence, if there existed a gravitational 
theory of cosmic-ray phenomena in the substratum, it would arise 
because a free particle which passes one observer with a small velo- 
city passes another observer with a velocity less than c but com- 
parable with c. However, we have found that, in whatever way the 
observers graduate their clocks, the velocity of a free particle is 
the same relative to each observer it passes, at the instant of passing. 
Consequently, there is no intrinsic acceleration of a free particle in the 
gravitational field of the ‘pure’ substratum and hence no gravitational 
explanation of cosmic-ray phenomena. It is still possible that large 
gravitational accelerations may occur in the real universe, owing to 
the presence of the nebular condensations which distinguish it from 
the ideal ‘pure’ substratum. 
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_ , Se 8 
[HE equation (scat ial ae = (1) 
occurs in a number of physical problems, such as that of the vibra- 
tion of membranes, when the boundary condition (at a supported 
edge) is w= 0; (2) 
and that of convection in a layer of fluid heated from below, the 
boundary condition then being 
~ = 0, (2a) 
Ov 
where év is an element perpendicular to the boundary. An exact 
solution of (1) has been obtained for one or other boundary condition 
in relation to several shapes of boundary, and for both when the 
boundary is a circle; but Rayleigh, in a paper concerned with the 
thermal problem,* remarked that solutions were unknown to him 
for boundaries having the form of a regular hexagon or of an equi- 
lateral triangle. 

Results obtained in an approximate (‘relaxation’) solution of (1) 
and (2a) for the hexagon suggested that an exact analytic solution 
does in fact exist; and, since on general grounds the mode w may 
be expected to be expressible in trigonometric functions, the fol- 
lowing results were obtained without difficulty. 

For a hexagon having sides defined by the straight lines 


V3 1 1 
(+ Sa)(etytn)(52-v8) = 0, (3) 


the required function is given by 
4nay 


9 9 
w= ju, 008 Sp (Var +y) +008 = (v32—y) +08: 3h (4) 


where w, is the amplitude of the vibration of the central point, and 
A in equation (1) is given by 


4nr 
A= —. 
3h 


* Phil. Mag. 32 (1916), 529-46; Collected Papers, 6 (1911-19), 432-46, 
No. 412. 
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This solution has some interesting properties. In the first place, 
since it has complete symmetry about the diagonals, it is a solution 
of the same problem for the equilateral triangular boundary. 
Secondly, the values of w on the lines 


3 1 v3,\/f/ 1.,% ‘ 
(v3) (2+ Wits a)(a — It a) == © (5) 


when n is odd are every wherbamortatgng, gnd are equal to dw). The 
position of these lines is indicated in Fig. 1, in which the dotted 
curve, approximating to a circle, represents the node in the case 
when 7 equals one. 


y 





























Fic. 1 


Finally, transferring the origin to a vertex of the original hexagon, 
for example the point (0,), we obtain 
9 1) 


ae Yt v8e—h) + 


9 
wW = —§Wp\COs 


(6) 


9 y 

2n7 ' 4n7 

+cos 3] (y—~v3x—h)-+-cos — 
) 


3h 


which gives w a constant value equal to 3w, round the sides of an 
equilateral triangle whose centroid is the origin, the central displace- 


ment again being wy. 
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The similar function 


9 9. 

w= a wo{sin = (v¥3a—y)-+sin se (v3a-+-y)-+sin a). (7) 
for which the value of A is the same as before and is again odd, 
satisfies (1) together with the boundary condition (2) on the hexagon 
(5), but it also vanishes on its diameters. It is therefore the solu- 
tion for an equilateral triangle subject to the boundary condition 
(2), i.e. it represents one possible mode of vibration for a triangular 
membrane having fixed edges. It does not appear possible to obtain 
in this manner the gravest vibration of the hexagonal membrane, 


and the approximate solution (obtained as mentioned above) gives 


no indication that any such solution exists. 








ON DIRICHLET’S SINGULAR INTEGRAL 
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1. Let f(t) ¢ L, over (—o, 00), let 
- : 1/p 
fly = ( [ or it (1<p<«) 
or |f |, = essential upper bound |f(¢)| (p = 0) 
—ao<f{<a 
respectively, let « be finite and positive throughout, and let D,f be 


Dirichlet’s operator 
sin a(t—s) 


Df = $16,0) =~ { so we. (1) 


Plainly Dirichlet’s integral exists in the ordinary Lebesgue sense and 
is uniformly bounded and uniformly convergent in —oo < s < o, 
when 1 < p < © since, by Hélder’s theorem, 


; re isin at , a 
oat <u] (i<r=2 <a), am 
p’ 


t p— 1 


Us . 
ron sin “(t—8) ay) di 
wT) t—s 


—U, 
=~, 


r \1/p os \ I/p 
<= : ( [ if ae) +( | ror a (U,, Uz, > 0). 
p’ U 


ny 
— 2 


Dirichlet’s integral is important not only in the theory of Fourier 
series, but also as a ‘resolution of the identity’, and as a matter of 
fact a great many properties of D,f were well known under certain 
assumptions on f(t) and not merely the assumption f(¢) ¢ L,. Some 
years ago E. Hille and J. D. Tamarkin* proved some further pro- 
perties of D,f for 1 < p < o, and they also gave an example which 


* Bull. American Math. Soc. 39 (1913), 768-74, referred to as H.-T. I, 
and Compositio Math. 1 (1934), 98-102 [H.-T. IT]. 
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showed that D,f does not necessarily belong’ to L,(—0o0,00) when 
p = 1. The main tool that they employ is Hilbert’s operator 


fle) = TH = 1 him [aes Se—) 2 fe >0). (1.3) 


TT «—0 


Hille* has also put forward some questions and conjectures. In this 
paper we shall deal with these and with some connected problems; 
and we shall discuss f(s,«) = D,f as a function both of a real and 
of a complex variable (§§ 2-4 and §§ 5-6 respectively), and we shall 
see that the latter case is adapted for elucidating some of the pro- 
blems on a purely real variable. As corollaries I shall give some 
theorems about integral functions and among them an extension of 
a well-known theorem of Paley and Wiener,+ based upon the repre- 
sentation of integral functions of a certain class by their proper 
Dirichlet integrals (Theorem 5, § 5). 


2. The case of a real variable 

The following theorem is valid: 

THEOREM 1. Let f(t) € L,( ,00) and let Df = f(s,x«) be the 
Dirichlet integral of f(t) 

(la) When p = 1, then D,f belongs to L, for allr >1 

(1b) When 1< p < ©, then D,f belongs to L, for all q > p, and 

(1b’) [Df lp < Co lf lp, where C,, depends on p only, 

(1b6") |D,f- fly —> 0 ( (~ > 00). 
(1b”) Let A, be the best possible value of C,, in (1b’), then A, = Ay 


[p’ =p toi, A,, 1s a monotone non-decreasing and continuous 


* E. Hille, Trans. — Math. Soc. 39 (1936), 131-53 [Hille I], and 
Proc. National Ac. of Sci. 24 (1938), 159-61 [Hille IT], and Annals of Math. 
10 (1939), 1-47 [Hille III]. A great part of the difficulties will be overcome 
in this paper by the equation (2.4), which has been known hitherto only for 
f(s) € Le, g(s) € Dg. 

+ R. E. A. C. Paley and N. Wiener, ‘Fourier transforms in the complex 
domain’: American Math. Soc. Coll. Pub. 19 (1934), 12-13. 

t H.-T. have proved (16) for 1< p< 2, gq=~p, (1b’) and (16”) for 
l p < 2, but their proof also holds for 2 < p < o, as they state in H.-T. II. 
Hille [I] has proved that A, =1 and has conjectured that A, > 1 for 
|< p < 2 (vide 1b”). This conjecture certainly is true when p is large and 
when p—1 is sufficiently small. Hille [I] has proved (lc) for 1< p < 2. 
From this case we shall deduce the general result, but there is no difficulty 


in starting from p = 2. 
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function of p (p > 2), and ap < A, < 4p/log2* for some a> 0 
(p > 2). 

(lc) When 1< p < @, then 

D,(Dgf ) = Dg(Dif) = Df [y = min(a, B)]. 

(ld) Le l<p<o,letl<r< oand rp < p+, let fe L, and 

g € L,; then . - 
| f(s, w)g(s) ds | f(s)g(s, x) ds, (2.4) 


(1d’) , f(s, «)g(s, x) ds S(s)g(s, «) ds. 

The proof of (1 6”) is left to the third section. Some other theorems, 
on the solutions of D,f = f, D,g = f, and D,f = 0, will be given in 
the fourth section. The proof of (16’) and (16”) was given by Hille 
and Tamarkin. That D,f belongs to L, for any q > p immediately 
follows from D,f¢ L, and from D, fe L,, (1.1). 

For proving (14) and (1d) we need the following 


Lemma 1. Let feLl,,geL, (l< p, r < 0), and let 


s—t 


f(s, «, U) -2 fe nF 1 


oo U 
Then | fe, a, U)g(s) ds = [ S(tg(t, «) dt. 
ae +u 
We have 

U U 0 ; 

[ sate, a) ae =~ f feo ae [ g(a) 209 as, 

-U -U ; 
and we can interchange the order of integrations in consequence of 
absolute convergence, using Hélder’s re mae ; hence (2.2) is true. 

Now let 1 < p< 0,1<r< wand let pr < p-++, that is to say, 

r<p' = p/(p—1). Since, by (15), g(t, x) belongs to L,,, the integral 
on the right-hand side of (2.2) converges for U +00; hence the 


— oo 


integral ” 
I(a, U;g) = f(s, «, U)g(s) ds 


* The corresponding results for Hilbert’s operator are due to E. C. Titch- 
marsh, Math. Zeits. 25 (1926), 321-47, Theorem C, and to M. Riesz, Math. 
Zeits. 27 (1928), 218-44 [Riesz I], and Acta Math. 49 (1927), 465-97 [Riesz IT]. 
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tends to a finite limit (U > oo) for any ge L,, hence f(s,«, U) con- 
verges weakly* to a function ¢(s) € L, [r’ = r/(r—1)], 


Het :q) > f $(s)g(s) ds (U-> 0), (2.3) 


and, since, by (1.2), f(s,«,U) converges uniformly to f(s,«), this 
function must be equivalent to ¢(s). So from (2.2) we get (2.4) for 
l<p<oao,l<r<o, prc pr. 

When p = 1, then r can be any finite number greater than unity; 
so can r’, hence the function f(s, «) = ¢(s) belongs to L,, for all finite 
r’ > 1, and by (1.1) plainly f(s, «) belongs to L,, also; hence (1a) is true. 

When 1 < p < ©, r = 1, then, by (1.1) and (1.2), 


[ {f(s, 0, U)—f(s, «)}g(s) ds 


a 


—-N 
( f+ f )tlste.0, C+ 106, a)B iat) as + 
- N 


- 
20 


+ [ f(s, 0, U)—fle, «)]Ig(s)| ds < « 


(V > Me), U > U,(N)); 


for f(s,a,U) is bounded uniformly in s and U, and f(s, «) in s, and 
also f(s,a,U) tends to f(s,«) uniformly in s. Hence, by (la), from 
(2.2) we get (2.4), and (1d) is true. 

Now (2.4) may be written in the form 


x 


f()Dg) dt (l< p,r < 0; pr < p+r). 


[ g(t)D,f dt = 


— 


Take r = min(p’, 2); then 1 < r < 2, and, as Hille proved, 
D, Dg : Dz Dg Dig (a, B > 0; y = min(a, B)); 


* The theorem on weak convergence is cited by 8S. Kaczmarcz and H. Stein- 
haus, Theorie der Orthogonalreihen (Warszawa-Lwow 1935). It follows from 
S. Banach’s Theorem 6, Théorie des opérations linéaires (Warszawa 1932), 80, 
and from the well-known ,theorem: When 1 <q< © and |z,(t)|, <A 
(n 1, 2,...), then the sequence {x,(t)} contains a subsequence which converges 
weakly to a function x(t)€ L,. In order to prove that f(s,a) = 4(s), take 
g(s) lforO0O<s<yory<s <0 (y= 0), g(s) = 0 otherwise, and make 


use of (2.3) and (1.2). 
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hence, in consequence of (2.4), 
| (D,, Dz f—D,f)g dt = | (Dz fD,.g—fD,g) dt 
- | (Dz D,g—D,g)f dt = | f()(D,g—D,g) dt = 0 


for all ge L,. 


7 


Thus D, Dgf—D,f = 90, and we have proved (1c). 
Taking y = 8 = a, we easily get (1d’) from (1c) and (1d). 

3. Proof of (1b”). Let 1< p< @ and g(t) € L,,; then in con- 
sequence of Hélder’s inequality from (2.4) and (16’) we have 


| f(s)g(s, a) ds < If(s, x) |» I9\p . A volt lolg ly 


7 8 
for all fc L,,; hence, by a well-known theorem, |g(s, «)|,, << A,|g|y, 


3 <A,, and so 


and consequently A,,<A,. By reciprocity A, 


A, = A,. Now log A, is a convex function of p (Riesz IT), 
— 2-2 


Per A,! “PAT? (l<p<r<q< 00). 


; 
Taking q > 2 and p = q’ = q/(q—1), we immediately get A, < A,; 
therefore A, does not decrease for r > 2, and plainly it is a con- 
tinuous function of 7. When 7'f denotes Hilbert’s operator and when 
p is an even positive integer,* then |7'f|, < (2p/log 2)|f|,, and sot 
| Daf |, < (4p/log 2)|f |p. By convexity, A, < 4p/log2 when p > 2. 
To prove A, > ap (p > 2) for some positive a, let 0 > p > 2, 
h=p'—l,a=1, 


1 (Q<t<h), 
and let ae | 
wei fo | 0 (otherwise). 
Then |f |, = h¥”’, 

h 

oe bo t) 
2 s+ 

0 a 

where (s) - h sins 


78 


A(s) = f(—s, «)—g(s) = ad cos u—sin u)(h+s—u). 


J. 


* Riesz I. There f(t) is real-valued only. 
+ Cf. H.-T. I, inequality (8). 
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' |u cos u—sin u 

We have |7A(s)| < h?max | ——— ; 
, 2 

s<ussth U | 


and therefore 
|A(s)| < Ah%(|s]+h) (|s] <1), |A(s)| < Ah?|s|-* (|s] > 1) 


for O0< h<}(l< p’ < $); hence: 


: 1 [ ; 
Ah*! h?P 


|A(s) |?’ ds < [ + [+f Aww +2" cas, 
; p’—1 p—tl 


. 


x - “oo 


h\sins 


| f(s, x) |p _— lf(—s, a) | »’ => me —Ah® pW 


8 |p’ 


. If (s,« Mp - S > (pryun’! sin 8 —A(hh)2lp 
| f(s) $)|p° m| 8 


p’ 
Now plainly 0 < c < (h*)¥”’ < C for 1 < p’ < 3, and 


oo (n+i)7 po 


in |p’ 
sms)" dt > 2-P'l2 > | |s|-”’ ds = 21-¥'2 p> 
8 
R=—O cnt ayer n=0 
1—p’93p’/2-1 
qi-P'2 . a 
a ts (1—3-*4.5-*,..) > — 
p—l p—l 
for some a > 0 when h = p’—1 is sufficiently small, since the series 


tends to $ for h->0. Therefore 


1/p’ : 1/p’ hi h|p’ 
A, = Ay > : ne —AC* a = : p—AC? > ap 
a7\p —1 Tp 
for sufficiently large values of p; but, since A, > A, = 1, we have 
A,, > ap for some a > 0 and for p > 2. This completes the proof. 
4. Let l<r< 2, let d(t) run through all functions belonging to 
L,.(—oo, 0), let IN” denote the limit in the mean of index 7’, and let 


u(x) = F(t) denote the lait as 


w(a) = (2 +f d(tje- dt (r’ = r/(r—1)), (4.1) 

x -N 
then F is continuous. Let 2 denote the set of these functions ¢(x) 
it is well known that £” is a sub-set of L,.(—oo,00) and different 
from DL, when 1< r < 2, and it is evident that Z™ contains all 


step functions. Let 


$(x) = Fit) (f(t) € B”); 
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it is well known* that this operation is not continuous when r’ > 2, 
and that N 
$(—z) = M(2n)- fu b(t)e*** dt. 
= -N 

The following theorem holds: 

THEoREM 2. The function f(t) belongs to L,,(— 
or D, f = 0 if and only if 

(a) f(t) has the form 

a 


f(t) = (22)-4 [ p(uje du (b(u) e B®), — (4.3) 
~ ae 


or fit) = am0(2n)-1( | - [aye du (W(u)e B®) (4,4) 
N->o ae - 
respectively, when 1< p < 2; 
(b) f(t) has the form 
f(t) = (2n)3 J d(uje du (p(u) € Ly), 


oy 


or S(t) = M(2n)- (J J - {)e (uje*™ du (d(u) € L,,) (4.6) 


No 


respectively, or if f(t) is the limit in mean of index p of a sequence of 
functions like those respectively, when 2 << p < «; 

THEOREM 2(c). Let f(t) € L,(—, 0); then D, f = f or Df = 0 if 
and only if the Fourier transform we) = e is null outside or inside 
the interval (—«, x) respectively. 

The proof of Theorem 2(a) has been outlined in substance by 
Hille (I); to prove that the conditions are sufficient, we have still 
to show only that the functions f(t) defined by (4.3) and (4.4) belong 
to L,. Now the sum of these functions equals F-' which plainly 
belongs to L,,, and the function f defined by (4.3) also belongs to 
L,, in consequence of (4.2). 

Proof of Theorem 2(b). Let 2<p<oo and f(t)¢ E™; hence 
f(t) = Fd (pe L,), f(t+s8) = F{d(x)e-**}. We have t“sint € L,, 
(27)-* (|z]) <a 


F(7—'t—sin at) = 0 (|x| > 
x a 
a 


hence f(s,a) = | fe+9™ x dt = (2m)-* | h(a)e-8% dx (4.7) 


—a@ 
* E. C. Titchmarsh, Proc. London Math. Soc. (2) 23 (1924), 279-89; H.-T. I. 
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by Parseval’s theorem. Therefore, from D, f = f and (4.7), it follows 
that f(¢) has the form (4.5), from D,f = 0 by the uniqueness theorem 
that d(x) must vanish for |z| < « and that f = F¢ has the form 
(4.6); evidently these conditions are also sufficient. Now let f(t) be 
an arbitrary function of L,, such that D,f = f, let {g,,(t)} (n = 0, 1...) 
be a sequence such that g, ¢ E™ and |g,—/f|, > 0 (n > 00), and let 
f, = Dyg,. Plainly D,g,, = 9,(8, «) also belongs to B® by (4.7), and, 
since D, f,, = DD gn) = DyGn =f, by (le), f, is of the form (4.5), 
and, in consequence of (1 5’), 


lf—Sale = IDA f—In) |p < CrlF—Inlp >0 (n ae 0). 
When D,f = 0, take g, ¢ E™ and |g,,—f |, > 0 (n > 00), and put 
Sn = In—DyG»; then f,, ¢ E™, D,f, = 9; hence f,, has the form (4.6), 


and the sequence {f,,} approximates to f, since 
lf—Saly _ lf—9.—D, f+D, In\p 
—_ lf- gal,+ |D,( (f— In lp (1+A,) \f- Inlp 
and |f—g,|,—> 0. The converse is quite simple. So Theorem 2 (5) 
is true. 
In Theorem 2 (c) less is proved than in 2 (a); it immediately follows 


from the equation 


a foe) 


= —- | 10 dt | etus—) dy = = | es" du | toe dt 


—a@ —@ 


ee. isu 
_ ay J oem du, 
the interchanging of the integrations being justified by absolute 
convergence. 

We next state the following theorems: 

THEOREM 3. Let f(t)e L, (l< p < 0); then D,f equals f if and 
only if f has the form f = D,g (ge L,); and Df = 0 if and only if 
f has the form f = g—Dyg (g € L,). 

THEOREM 3(b). A necessary and sufficient condition for the equation 
Dg = f to have a solution g € L, (1 < p < ©) is that the function f(t) 
belongs to L,, and is a solution of D,f = f. 

The proofs rest on Theorem 1 (c). When p = 1, in Theorem 3 (5) 
the condition f¢ L, is not necessary, but is one of the sufficient 
conditions; possibly it can be replaced by a certain weaker condition, 
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which plainly is necessary, namely that f(t) belongs to L, for all 
> 2 


5. The case of a complex variable 

Let C denote an arbitrary straight line in the complex plane 
extending to infinity in either direction and making an angle ¢ with 
the positive real axis; when ¢ = 0, then C either is parallel to the 
real axis or coincides with it. We again take 0 < a < o. The fol- 
lowing theorems are valid: 

THEOREM 4. Let —l<r <1 and f(z)(1+|z|)’ ¢ L, along C; let 
B = ae-*?, Then 


fls,B) = 1 tz) (2) — ) dk 


is an integral function of s and 
Lf(s, B)| < ’ 1+ |R(se-) |} exp|aB(se-*)]. 
THEOREM 4(a). Let f(z) belong to L, (1< p < @) along C, then 
g p I g 
f(s, B) (8B = ae-'*) is an soot function of s and 
Lf(s, B)| < A exp|a3(se~4)], 
THEOREM 5. Let 
(i) f(z) be an — function, 


(ii) a Him max log fl M be finite, 
r 


(iii) Ste) z) Fy {)-* baa to L,, along C for some p (0< p< 1), and 
let argB = —¢. Then, for Be’? > x, 


= | fe) a $) de = fis). 


This identity is not true when Be'? is less than ap. 


THEOREM 5(a). The preceding theorem also holds when we replace 
(ili) by the condition (iii’) f(z) € L,, along C for some finite positive p. 


THEOREM 6. A necessary and sufficient condition for the function 
S(tjeL,(—~,0) (l< p<) to be a solution of the equation 
D, f = f(t) is the existence of an integral function f(z) which satisfies 
the condition |f(z)| < A,e+©| for all « > 0 and which equals f(t) 
almost everywhere on the real axis. 
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THEOREM 7. Let 
(i) f(s) be an integral function of s, 
(ii) | f(s)|< A, e+! for all « > 0, 
(iii) —1<r< 1* and f(s)(1+|s|)’ belong to L,, along C (0 <p <1). 
Then |f(s)| < A{1+|R(se-#)]} “ exp|a3(se-#)|. 


= 


THEOREM 7(a). When (iii) is replaced by the condition f(s) € L, 
(0 < p < o) along C, then |f(s)| < A exp|aS(se~?)]. 


THEOREM 8. Let (i) f(s) be an integral function, (ii) | f(s)| < A, e&*+©*! 
for some « > 0 and all « > 0, (iii) f(a) € L,(—, ©) for some finite 
p > 2. Then either the function f(c) is the Fourier transform of a 
function $(t) € L,,(—0, 0) which is null almost everywhere for |t| > «, 
or else, if f(a) is not a Fourier transform, then f(c) is at any rate the 
limit in mean of index p of functions f,(c) = Fd, (n =1, 2,...) of 
this kind. 

The conditions (i), (ii), (iii) are also necessary. 


6. When proving these theorems we take C as the real axis, the 

general case being easily deduced from this one by the substitutions 
zo t+wel?, s = z)+ue'?, w real along C. 

Proof of Theorem 4. We have ¢ = 0, 8 = a; the integrand in (1) is 
an integral function of s = o+i7, and the integral converges uni- 
formly in any circle |s| < R; therefore f(s) is an integral function. 
Now we put 


sin a(t—s) 


+1, 
J a= htt 


<r <0. Since |t—s| >1 in J, and J,, we have 


1+ |¢| 


; t—s| > |t—s|-" > (2+]e|)"(1+|é])~, 
2+ |o| gia 


|t—s| > |t—o| > 


so then in consequence of the inequality |sin a(¢—s)| < e™"! and of 


the above hypothesis 
[I+ lel < A(2+ lol)Fet! < 2-7A(1+ Jol) rer" 
* When r > 1, then evidently the condition on f(s)(1-+- |s|)~1 is also satisfied. 


When r is any positive integer we can get a slightly better result than that 


stated in Theorem 7. 
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In J, we have |t—o| <1, |é] < |o|+1, and so, for -—l<r< 0, 
1 < (2+|o|)(1+]é))* 
1 < (2+ |o|)7(1+ |él)” < 27(1+ lo|) "(1+ [a] 


o+1 ro) 


[Ig] < Ae! [ | f()| dt < A2e*"\(1+ Jol) | [f(Q|(1+ [ely de 
o—1 — 
because |sin a(t—s)/(t—s)| < Ae*!; and 
If(e,0)| < A(1+|ol)-rex" 
now follows for —l1<r<0. When 0<r <1 we get the same 


result by making use of the inequalities 
lt—s| > |t—o] > 271+ for + Ie)" (|t—o] > 1 
L< (2+(t\+ ol)" < 21+ [e+ ol)" (t—-o] <1 
To prove Theorem 4(a) we make use of Hélder’s inequality 


M sin a(o+i7r—t,) ., | (sinh2a7r-+-sin?a(ao—t)) 4 
a|f(8,~)| < lf |p “| =Iflp . 
o+ir—t » (o—t)?+7? 9 


me sinh o sin a(o—t) 
IF l| nel + =e = lf lof 4+}. 


{( —t)?+77}4 p’ o—t p’ 
Now J, is a constant which depends on « only; when p > 1, 
o+1 , : o--1 © ;, 
Bs sin p y sinh ar?’ i 
W< | SS) @44 4 + et dt < Aer'atr, 
: J o—t 
go=—1 — © o+1 
so that J, < Aexp|at|, and this inequality also holds when p = 1. 
Therefore | f(s, «)| < Aexp|a7|, which completes the proof. 
Proof of Theorem 5(a). In consequence of the 
Lemma 2.* W baler is an integral function and ey )| < Ae! and 
f(t) € L,(—, ©) for some finite p > 0, then f(t) > 0 (t+ +00) and 
SiH «€ oe 00) for anyq > p 
we have f(t) > 0 (t-> +00) and, by Hélder’s inequality, 


FHF é|)-? € L,(—oo, 00). 
Hence f(s, «) is an integral function, and we may take s real without 


* M. Plancherel and G. Pélya, Comm. Math. Helvetici, 9 (1936/7), 224-48, 
and 10 (1937/8), 110-63; vide I, §§ 15-17, and IT, § 30. 
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loss of generality. First, we shall show that, for |y| > a)+8 (8 > 0), 
the function* 


N 
dy(y) = (27)-4 | flo)e-*2" do 
-—-N 


converges to zero uniformly when WN tends to infinity. 
Let « = 48 and £ = a9+2c, and let G,(z) = f(zje* (2 = u+iv), 
then 
IG¢(z)| < | f(z)|ebi2! < A, dvort+eiel = A,e, 
|Ge(u)| = |f(u)| +0 when w>+o0, 


|Ge( iv)| = [f(a lege <f A, ror ellvl+ge +>0O when v—>-—oo. 


Therefore, by applying well-known theorems of Phragmén-Lindelof ¢ 
to Gz in both quadrants of the lower half-plane, we see that the 
function p(r) = max|G@;(re'*)| (—7 < ¢ < 0) is bounded and that 
it tends to zero for r > oo. Now 
N N 
| f(c)e-*v do = f & o)ei*t-—v) dg = iN fee Gi( Nett )etl-wne*eid dd, 
—N —N 
and so for y > ag +8 = ap +3e = (+e 
0 
z)*|by(y)| < w(N)N | eNw-dsing dd 
im ia 
= 2u(N) | eN(C-y)sin d VV dd x 2u(N) | e-Ne2dlaNV dd < e-lnpy(N), 
0 0 
so then dy(y) > 0 (N + o) uniformly for y > a +5; the correspond- 
ing result holds for y < —a )—8. 
Now, in consequence of absolute convergence, 
B N B 
Qn) [ py(y)eiu dy = (2x)? [ fit) dt [ etsv-itv dy 
+ yy 4 


= 1 { no2ee— dt. 


t—s 
—-N 
* Compare the similar investigation by Plancherel and Polya I, § 19. 
+ E. C. Titchmarsh, The Theory of Functions (Oxford, 1932), (5.61) and 
(5.63). 














78 H. KOBER 


The right-hand term tends to f(s, 8) for N > 0. Let B > a +8; then 
the left-hand side equals 


7 uy 6 es Oe— ry B 
(2m) [ dy(ye dy tml f + f ldn(yeiu dy 
. | j : J 
—a%—8 —B Pe 


= it Lh+ 4; 
and J, 0, I, > 0 (N +0); and so the limit of J, (VN > 0) exists 


also and equals f(s, 8), 
+8 


(22r)-* lim [ dy(y)e*” dy = f(s, B) 


N+o _,*_5 
for any B > a +8; but here the left-hand side does not depend on 
B, and f(s, 8) converges to f(s) for B > 0, since f(s)(1+-|s|)-? belongs 
to L,(—0oo,00) and f(s) is continuous.* Hence f(s, 8) = f(s) for any 
B > a»; now f(s, 8) is a continuous function of 8 when s is fixed, and 
therefore also f(s, a9) = f(s). Incidentally, we see that f(s) must vanish 
identically when a, = 0.T 

If f(s, 8) were equal to f(s) for some positive 8 < a», then, taking 
s complex, by Theorem 4(a) we should get the contradiction 

|f(s)| = | f(s, B)| < Ae”* for y= }(at+B) (0<B<y < o%). 
The case 8 <0 is now evident, since f(s,8) = —f(s,—f). So 
Theorem 5 (a) is proved. 

Proof of Theorem 5. Let z, be a zero of f(z); if there were no 
zero, in consequence of Hadamard’s factorization theorem, f(z) would 
have the form bexp(az); but (1+ ]|¢|)-texp(at) does not belong to 
L,,(—«,«) for 0 < p <1, and so we should get the contradiction 
b = 0,f = 0. Now &(z) = (z—z,) f(z) is an integral function, plainly 
|ys(z)| < A, exp|(a )+e)z| for any « > 0, and, by (iii) %(t) € L,(—%, %) 
(0 < p <1); hence, by Theorem 5 (a), we have 


sin B(t—s) 


dt, 
t—s 


wa) == | oe 
and so ies 


f(s) = $(s)(s—z.) = : | fey eB) ay 


@ 
t—s 7 


[ b(t)sin B(t—s) dt 


oi (5.1) 


* E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals (Oxford 
1937), Theorem 12. 
+ Cf. V. Ganapathy Iyer, Quart. J. of Math. (Oxford), 7 (1936), 294-9, 


corollary 2. 
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for any B > a. In consequence of Lemma 2 we have (t) € LZ, and 


b 


[ b(t) dt | cos B(t—s) dB = f(s) 


b oo 

api | ib(t)cos B(t—s) dt = : 
F 7 E & ak 
for b > a». Now the inner integral on the left-hand side is a con- 
tinuous function of 8 and must therefore vanish for 8 > a, identically 
in 8; consequently the integrals 


[ sb(t)sin Bt dt, | ih(t)cos Bt dt, | i (t)sin B(t—s) dt 


vanish* for B >a». So from (5.1) we get f(s) = f(s, 8) for B > a, 
and the theorem is proved. 

When 8 < a, we use the same argument as in the proof of 5(a). 
Of course, the theorem is true also when the condition (iii) is replaced 
by the stronger condition 

fOHU+I|é) € L,(—0,0) (0O<p<cl;—l<r<l). 

From the Theorems 4(a) and 5(a) we immediately get another 
proof of the assertion Dz D,f = D,f made in the Theorem 1 (c) for 
the case 8 > «a. Theorem 6 follows from the Theorems 4 (a) and 5 (a) 
at once; so does Theorem 7 from 4 and 5. 

The Theorems 2 and 6 give two different answers to the same 
problem, that of the solutions of D,f =f. By comparing them we 
get another result which, for p = 2, is a well-known theorem due to 
Paley and Wiener, while, for 1 << p < 2, it is due to Plancherel and 
Pélya.t For p > 2 the result has not previously been given; so for 
2 < p < o this Theorem 8 follows. 

The part played by the line C in some of the preceding theorems 
is elucidated by the following: 

Lemma 3. Let f(z) be an integral function and |f(z)| < Ae®*|, let 
m be an integer and f(z)(1+-|z|)™e L, (0 < p< @) along C. Then 
f(z) A+ |z|)™ € L, for any q > p along every line C’ parallel to C. If 
there is another line C, which is not em to C and along which 
f(z)\(1+|z])™ belongs to L, (0 <r < a), then f(z) must vanish identi- 
callyt when m > —l, while it vanishes identically or reduces to a poly- 
nomial of degree less than |m|—l when m < —l, wherel = max(1/p, 1/r). 


e 


0 


This is Paley-Wiener’s theorem for p = 1. 
+ Plancherel and Pélya I, §§6 and 15. The result follows from that of 
Paley-Wiener by Lemma 2. 
t Cf. V. Ganapathy Iyer, loc. cit., corollary 3. 
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The proof of the first assertion rests on Lemma 2 and on the 


theorem: 
Let f(z) be an integral function but not a polynomial, and let 
| f(z)| < A, exp|(a+e)z| for all « > 0, let m be an integer and 


fl(1+|t))™ € L,(—20,00) (0<p<o), 
let 21, Za,.--5 Zim) be arbitrary numbers or zeros of f(z) when m > 0 or 
m < 0 respectively, and let g(z) = [J (e—z,) (k = 1, 2,..., |[m]), TT =1 
(m = 0); then 
M co 
[I flcbiyg(a-iyymm|? dex < ex’ [ | f(x)g(a) mr dex. 
—w a 
This theorem was proved by Plancherel and Pélya* in the case 
m = 0; the general case easily follows by an argument similar to 
that used in the first part of the proof of Theorem 5. Under the 
hypothesis of Lemma 3 plainly f(z) has an infinite number of zeros 
when m > —1/p, while f(z) reduces to a polynomial of degree less 
than |m|—1/p when m < —1/p and f(z) has a finite number of zeros 
only. The second part of Lemma 3 rests on Lemma 2 and on the 
theorems of Phragmén-Lindelof cited above. 


* Plancherel and Pélya II, § 27. 








